KSU   PILE   STANDARDIZATION   AND  STUDY    OF 
SLOWING   DOWN   AND   DIFFUSION   MODELS 


by 


LARRY    RAY    FOULKE 


B.   S.,    Kansas  State  University,      1960 


A   THESIS 


submitted  in  partial  fulfillment  of  the 


requirements  for  the  degree 


MASTER    OF   SCIENCE 


Department  of  Nuclear  Engineering 


KANSAS   STATE    UNIVERSITY 
Manhattan,     Kansas 


1961 


LD 

Ti  t 
I  ft/ 

£  .     i_  TABLE   OF   CONTENTS 

INTRODUCTION 1 

THEORY  4 

Thermal  Neutron  Diffusion 4 

Neutron  Slowing  Down 6 

Solutions  for  Particular  Conditions 8 

EXPERIMENTAL  FACILITIES    16 

General  Pile  Description 16 

Standard  Pile  Configuration 19 

Finite  Pile  Configuration 19 

Foil  Counting  Facilities 22 

Neutron  Probe  Facilities 27 

Traversing  Mechanism 29 

EXPERIMENTAL   PROCEDURE 34 

Foil  Measurements 34 

Neutron  Probe  Measurements  in  Standard  Pile 36 

Neutron  Probe  Measurements  in  Finite  Pile 38 

DATA   PRESENTATION  AND  ANALYSIS 38 

Pile  Standardization  Data  Treatment 38 

Empirical  Fit  of  Indium  Resonance  Data 44 

DISCUSSION   AND  RESULTS 50 

Standardization  Constants 50 

Confidence  Limits  on  Flux  Assignment 57 


TABLE    OF   CONTENTS    (concl. ) 


Intel-calibration  with  Argonne  Standard  Pile  and 

Hanford  Standard  Pile 65 

Theoretical  Standardization  of  BF3  Neutron  Probe 69 

Scintillation  Probe  Standardization 70 

The  Use  of  Gaussian  Ranges  and  the  Fermi  Age  to 

Describe  Slowing  Down 71 

Study  of  Source  Boundary  Condition  in  Finite  Pile 73 

Suggestions  for  Further  Work 89 

ACKNOWLEDGMENT 91 

LITERA  i'URE   CITED 92 

APPENDICES 96 

APPENDIX  A:  Solutions  to  the  Fermi  Age  Equation  and  the  Thermal 

Diffusion  Equation  for  the  Standard  Pile  Analysis.  ...         97 

APPENDIX  B:  Finite  Pile  Solution  to  the  Thermal  Diffusion 

Equation  Using  the  Distributed  Thermal  Source 

Found  from  the  Solution  to  the  Fermi  Age  Equation...     102 

APPENDIX  C:  Finite  Pile  Solution  to  the  Thermal  Diffusion  Equa- 
tion Using  a  Point  Source  of  Thermal  Neutrons 106 

APPENDIX  D:  Description  and  Explanation  of  the  IBM-650  Compu- 
ter Program  Used  to  Calculate  the  Flux  in  the 
Standard  Pile 109 

APPENDIX  E:  Description  and  Explanation  of  the  IBM-650  Compu- 
ter Program  Used  to  Calculate  the  Thermal  Flux 
from  a  Fast  Point  Neutron  Source  in  a  Finite  Pile.  .  .       117 

APPENDIX  F:  Description  and  Explanation  of  the  IBM-650  Compu- 
ter Program  Used  to  Calculate  the  Thermal  Flux 
from  a  Thermal  Point  Neutron  Source  in  a  Finite  Pile.  127 

APPENDIX  G:  Description  and  Explanation  of  the  IBM-650  Compu- 
ter Program  Used  for  Empirically  Fitting  Experi- 
mental Data  with  the  Sum  of  Several  Terms  of 
Gaussian  Form 137 

APPENDIX  H:  Considerations  of  Perturbations  and  Corrections 
Encountered  in  Experiments  Using  Finite  Foils 
and  Sources 146 


INTRODUCTION 

In  investigations  concerned  with  neutrons  and  their  behavior,   there  is 
often  a  need  for  accurate  measurement  of  a  laboratory  standard  of  abso- 
lute thermal  neutron  flux.     In  1941,    Anderson  and  Fermi  (3)  outlined  a 
method  of  determining  an  absolute  thermal  neutron  flux  in  a  laboratory 
pile,   based  on  theoretical  and  semi -empirical  considerations.     Seren  (35) 
seems  to  have  been  the  first  to  do  the  experimental  work  in  his  calibration 
of  the  Argonne  Standard  Pile  in  1943.     This  pile  was  dismantled  in  1946 
and  Redman  (32)  recalibrated  what  is  now  the  Argonne  Standard  Pile.     In 
195  I,    Hanford  Laboratories  (11)  published  a  report  on  the  standardization 
of  indium  foils  and  gold  foils  in  the  Hanford  Standard  Pile.     It  is  reason- 
able to  think  that  other  piles  have  been  standardized  as  this  method  of  pile 
standardization  has  been  outlined  in  many  references  (11,    16,    20,    30). 

A  standardized  pile  is  a  large  stack  of  graphite  in  which  a  neutron 
source  of  known  strength  has  been  placed,    and  an  absolute  thermal  neutron 
flux  has  been  defined.     The  spatial  distribution  of  the  slowing  down  density, 
q,    at  indium  resonance  energy  can  be  measured  with  cadmium -cove  red 
indium  foils.     Assuming  no  leakage  of  epithermal  neutrons  in  the    z 
direction  of  the  assumed  infinite  column  and  no  absorption  of  epithermal 
neutrons,    absolute  values  can  be  assigned  to  this  measured  distribution  of 
the  slowing  down  density  since  the  volume  integral  of  this  slowing  down 
density  must  equal  the  known  neutron  emission  rate  of  the  source,   Q0. 
Empirically  fitting  the  absolute  slowing  down  density  leads  to  an  analytical 
expression  which  may  be  evaluated  at  thermal  energy,    and  which  then  can 


be  used  as  a  source  term  in  the  one  velocity  thermal  diffusion  equation. 
The  solution  to  this  equation  is  then  taken  to  define  the  absolute  flux.     It 
is  customary  to  use  the  neutron  velocity  of  2,  200  meters  per  second 
(m/sec)  in  the  calculation  of  the  flux  in  a  standard  pile  (20). 

The  confidence  limits  of  such  an  absolute  flux  assignment  are  deter- 
mined primarily  by  the  accuracy  of  the  source  strength  and  the  diffusion 
coefficient.     It  is  reported  by  Hughes  (20)  that  this  method  of  flux  assign- 
ment is  within  the  error  in  the  knowledge  of  the  absolute  source  strength. 
Previous  investigators  (11,    35)  seemed  to  have  no  question  about  un- 
certainty of  the  diffusion  coefficient.     However,    recent  measurements  of 
the  transport  mean  free  path  and  the  diffusion  coefficient  in  graphite 
introduce  additional  uncertainty  about  assignments  of  absolute  flux  (6  ,  29) 

The  National  Bureau  of  Standards  (7)  has  made  provision  for  pro- 
ducing a  constant  thermal  neutron  flux  of  4,  000  neutrons  per  square 
centimeter  per  second  (n/cmz-sec)  in  a  particular  moderating  geometry. 
This  calibration  has  been  made  to  an  accuracy  of  ±  2  per  cent  by  the  abso- 
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lute  counting  of    a    particles  from  the  B      (n,a)    Li  (7)  reaction.     It  has 
been  verified  to  within  the  limits  of  experimental  accuracy  by  a  compari- 
son with  an  independent  calibration  by  Oak  Ridge  National  Laboratories. 
The  National  Bureau  of  Standards  invites  the  exchange  of  activated  gold 
foils  between  it  and  other  national  laboratories  for  the  purpose  of  inter- 
comparison  of  absolute  flux.     This  has  not  been  done  at  Kansas  State 
University;     however  intercomparisons  have  been  completed  with  Hanford 
Laboratories,    and  with  Argonne  National  Laboratory. 

Glasstone  and  Edlund  (15)  have  assumed  that  for  the  experimental 
determination  of  the  diffusion  length,    extraneous  sources  to  the  thermal 


diffusion  equation  emit  only  thermal  neutrons.     This  development  appears 
to  be  adequate  for  describing  what  happens  to  the  thermal  flux  near  a 
concentrated  source.      The  method  of  assuming  that  the  source  emits 
epithermal  neutrons  of  age  zero,    and  allowing  their  distribution  in  space 
and  energy  to  be  determined  by  the  Fermi  age  equation,    leads  to  the 
more  accurate  expression  of  a  distributed  thermal  source  for  use  in  the 
thermal  diffusion  equation. 

Some  work  with  this  coupling  of  the  Fermi  age  equation  and  the  thermal 
diffusion  equation  in  addition  to  pile  standardization  has  been  done.     A 
group  at  Hanford  (10)  used  this  method  in  an  investigation  of  the  diffusion 
length.     They  concluded  that  generally  this  fast  source  theory  was  pref- 
erable,  but  only  in  the  case  of  a  single  source  did  the  diffusion  length 
obtained  from  thermal  source  theory  differ  appreciably  from  that  ob- 
tained from  fast  source  theory.     A  paper  by  Anselone  (4)  at  Hanford  in 
1954  outlines  the  use  of  the  "method  of  images"  in  determining  theoret- 
ically the  distribution  of  thermal  neutrons  from  fast  sources  in  exponen- 
tial piles.     A  publication  by  Wallace  and  LeCaine  (42)  has  an  extensive 
compilation  of  these  coupled  solutions  for  various  geometries. 

Previous  work  at  Kansas  State  University  leading  to  this  work  has 
been  the  determination  of  the  Fermi  age  both  from  theory  and  experiment 
by  Steichen  (3  9)  and  the  determination  oi  diffusion  length  in  the  KSU  pile 
by  Kaiser  (22). 


THEORY 
Thermal  Neutron  Diffusion 

Extensive  treatment  of  the  theory  of  thermal  neutron  diffusion  can 
be  found  in  many  references  (2,    15,    24,   43).     It  will  be  considered  suf- 
ficient for  this  work  to  outline  the  basic  theory  of  neutron  diffusion  and 
point  out  the  assumptions  and  approximations  involved. 

The  rate  of  change  of  the  number  of  thermal  neutrons  in  an  elemental 
volume  equals  the  number  which  flow  in  through  the  boundary  of  the  ele- 
mental volume  plus  the  number  produced  by  sources    minus  the  number 
captured  per  unit  time.     The  equation  expressing  this  principle  of  con- 
tinuity at  a  steady  state  condition  is 

-  V-T    +    S(7)  -    %n(?)  v  =  0,  (1) 

where    n    is  the  density  of  the  thermal  neutrons,     v    is  their 
representative  velocity.     £      is  the  absorption  cross  section  of 
the  medium,   and    S    is  the  number  of  thermal  neutrons  produced 
per  unit  time  per  unit  volume.     The  net  number  of  neutrons,     j    , 
flowing  in  unit  time  through  a  unit  area  normal  to  the  direction  of 
flow  is  the  neutron  current  given  by  Fick's  law  of  diffusion,  the 
characteristic  assumption  of  elementary  diffusion  theory 

f   =    -  DVnv.  (2) 

In  this,     D    is  the  diffusion  coefficient  for  the  monoenergetic  flux, 
having  the  dimensions  of  length. 


Using  Eq.  2  in  Eq.  1  gives  the  thermal  diffusion  equation,  Eq.  3. 
Since  for  monoenergetic  neutrons  the  neutron  flux,  (j),  is  equal  to  nv, 
this  equation  will  be  written 

D  V2  >|>  -Sa<t>  +  S  =  0  (3) 

where  an  isotopic  medium  has  been  assumed. 

The  conditions  for  the  justification  of  diffusion  theory  and  the 
diffusion  equation  include  the  following  considerations: 

(a)  Collisions  between  neutrons  do  not  occur,   making  the  pro- 
bability of  collisions  proportional  only  to    <)>.     This  results  in  linear 
equations  describing  the  neutron  diffusion.     Weinberg  and  Wigner 
(43)  point  this  out  as  the  most  important  mathematical  property  of 

the  diffusion  equation.     This  makes  superposition  of  solutions  possible. 

(b)  The  neutron  scattering  is  spherically  symmetric  or  isotropic. 
The  isotropic  distribution  allows  the  velocity  vector  to  be  treated  as 

a  scalar  in  the  equation  of  continuity.     This  assumption  places  a 
limitation  on  the  diffusion  equation  near  sources,   boundaries,    and 
absorbers  since  near  them,   the  scattering  is  not  isotropic. 

(c)  The  nuclei  of  the  isotropic,   homogeneous  moderator  are  at 
rest  and  of  such  high  mass  number  with  respect  to  a  neutron  that  they 
do  not  recoil  appreciably  under  the  impact  of  neutrons  (43). 

(d)  The  neutrons  are  monoenergetic  for  which  scattering  occurs 
without  energy  loss.  This  allows  the  diffusion  equation  to  be  energy 
independent. 


(e)  The  fractional  change  in  the  neutron  flux  and  the  mean  free 
path  over  a  distance  of  two  or  three  mean-free  paths  is  small  (15). 

It  can  be  shown  on  the  basis  of  transport  theory  that  Fick's  law  is  not 
valid  if  this  condition  does  not  hold  (15).     Thus  in  regions  near  a  con- 
centrated source,   absorber  or  boundary  where  the  flux  changes  rapidly 
or  the  mean-free  path  is  discontinuous,   the  neutron  current  cannot  be 
represented  as  being  proportional  to  the  gradient  of  the  flux. 

(f)  The  diffusion  medium  is  a  weak  absorber  so  that  the  scatter- 
ing cross  section  can  represent  the  total  cross  section,   and  therefore 
E    /  S      «  I.      This  is  tied  in  with  condition  (e)  in  that  if  the  absorp- 
tion cross  section  is  appreciable,    the  neutron  density  may  vary  sub- 
stantially within  a  mean-free  path. 

(g)  The  diffusion  coefficient  is  independent  of  position.     This  allows 
D    to  be  factored  out  of  the  term  for  divergence  of  the  current  in  Eq.    3 
leaving    DV2<(>    for    V-    DVf 

(h)  The  general  dimensions  of  the  moderator  are  much  larger 
than  characteristic  lengths  such  as  the  diffusion  length  in  diffusion 
theory  or  the  slowing  down  length  in  slowing  down  theory. 

Neutron  Slowing  Down 

The  problem  of  neutron  slowing  down  must  be  considered  in  order  to 
determine  the  effective  source  to  the  thermal  diffusion  equation.      The 
accurate  theory  here  is  quite  complicated  (43).     For  this  treatment  it 
will  be  sufficiently  accurate  to  consider  the  slowing  down  process  to  be 
described  by  an  approximate,    continuous   slowing  down  model  given  by 


Fermi.  It  is  based  on  the  assumption  of  a  functional  relationship  be- 
tween a  neutron's  energy  and  its  age,  i.  e.  ,  the  time  elapsed  since  its 
birth  at  a  primary  source. 

The  approximation  is  that  for  large  neutron  populations  the  neutrons 
lose  energy  in  a  continuous  fashion  rather  than  in  discontinuous,    dis- 
crete jumps  as  with  a  single  neutron.     This  approximation  cannot  be 
valid  in  the  event  that  a  neutron  may  lose  its  entire  energy  in  a  single 
collision  which  is  the  case  with  moderators  of  low  mass  number, 

The  development  of  the  Fermi  age  equation  which  describes  the 
neutron  slowing  down  process  is  given  in  considerable  detail  in  a  num- 
ber of  works  (2,    15,    43).     Assuming  no  capture  of  neutrons  during 
slowing  down,   the  Fermi  age  equation  is 

V'q(?,    r)=      3^f'   T>     ■  (4) 

A  new  variable,   the  Fermi  age  has  been  introduced  and  is  defined  by 

T(Eo,   E)     =     [E°      °(E)dE_    .  (5) 

*         £  Ss  (E)  E 

In  the  above  expression,   the  slowing  down  density    q(r  ,    r  )  can 

be  defined  as  the  number  of  neutrons/cm3  -   sec  which  slow  down 

past  a  given  energy  E,    corresponding  to  the  age     t,     at  a  given 

point     r.       E0  is  the  energy  of  the  source  neutrons,     E     is  the  energy 

to  which    t    is  measured  and     £    is  the  average  logarithmic  energy 

decrement  per  collision.     D  is  the  diffusion  coefficient,     and     2      is 

s 

the  macroscopic  scattering  cross  section. 


Solutions  for  Particular  Conditions 

The  solution  to  the  Fermi  age  equation  for  a  point  neutron  source 
emitting  Q0  monoenergetic  neutrons  per  second  in  a  finite  parallele- 
piped of  sides  2a,    2b,    and  2c  is  derived  in  Appendix  B,    and  given  here 


q(x,y,z,r)=£   J!      Y    QoFm(x)Fm(x0)Gn{y)Gn(yo)Ho(Z)Ho(z0)e-'*m+^tV'' 

(6) 


m  =  l    n  =  l    0  =  1 


where 


F     (x)    =    cos  a x  for    m    odd 

m  4S 


F     (x)    =    sin  a x  for    m    even 

nrr    '  I-  m 

Va 


F     (x„)    =    cos  a    x0  for    m    odd 

nv  /—  m 

va 


F     (x0)    =    sin  a    x0  for    m    even 

m  r~  m 

Kja. 


m  ir 
m  2  a 


G_(y)i    Q _(yo).    H   (z0)    are  defined  in  the  same  manner,   with  the 
dimensions    b    and    c    corresponding  to  directions    y    and    z, 
respectively.     This  solution  will  be  referred  to  as  the  slowing  down 
density  in  a  finite  pile  from  a  fast  source.     It  will  be  evaluated  at  the 
age  corresponding  to  thermal  energy  and  used  as  the  distributed 
thermal  neutron  source  in  the  thermal  diffusion  equation. 


(V) 


The  one -dimensional  solution  to  the  Fermi  age  equation  for  a  point 
source  of  monoenergetic  neutrons  at  the  origin  of  an  infinite  column  of 
rectangular  sides    a    and    b    is  derived  in  Appendix  A  and  given  here  as 

q(,,T).        4Q°     Y    Y     e-MiVa'  +  Wb^-zV^  (8) 

ab^  jti  k=l 

In  this,     z    represents  the  distance  from  the  source  along  the  central  Z 
axis  of  the  infinite  column,   and    j    and    k    are  summed  over  odd  integers. 

The  "7  T  is  often  called  the  slowing  down  length  or  characteristic 
length.  The  quantity  2\l  t  will  be  referred  to  as  the  Gaussian  range 
and  denoted  by 

r  =  2nTT    .  (9) 

Defining 

/r2 


y  y  e-4Mj2A2  +  k7b2). 


abWr    . 

j=l      k=l 

Eq.    8  becomes 

q(z,   T)=QoCe-z2/rJ. 

A  quantity  which  is  proportional  to  the  slowing  down  density 
may  be  measured  with  some  accuracy  with  cadmium -covered  indium 
foils.     There  are  some  correction  factors  that  must  be  applied  to 
this  measured  distribution  which  will  be  introduced  in  the  analytical 
part  of  this  work.     The  cadmium  covering  filters  out  most  of  the 
neutrons  with  energies  of  0.4  ev  and  below  so  that  the  activation  of 
the  indium  foil  is  due  primarily  to  neutrons  with  an  energy  of  1.44  ev, 


(10) 


(11) 
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corresponding  to  the  energy  of  the  resonaiice  cross  section  of  indium. 

From  the  solution  to  the  Fermi  age  equation  in  the  infinite  column, 
Eq.    1 1  ,    it  can  be  seen  that  the  natural  logarithm  of  the  activity  of 
cadmium-covered  indium  foils  plotted  against    z2    should  yield  a  straight 
line  with  a  slope  of  -  l/r2.       However,   the  neutrons  from  a  Pu-Be  source- 
are  not  monoenergetic  but  have  a  wide  energy  distribution  as  shown 
in  Fig.    1.     Hence,   it  is  not  unreasonable  to  find  a  curved  line  for  this 
plot  in  Fig.    14.     This  curve  might  be  considered  to  represent  the  super- 
position of  the  solutions  for  an  infinite  number  of  sources,    each  having 
a  different  monoenergy  and  thus  a  different  Gaussian  range.     Thus 
Eq.    1 1    should  be  written 

,2/r2 


v  - z  /ri 

q(z,r)  =  Q0  £      *j    C.  8 


(12) 


i  =  l 

where    f.     represents  the  fraction  of  the    Q0    neutrons  emitted  per 
second  that  are  in  the  i       energy  group    and    r.     represents  the 
Gaussian  range  from  the  energy  level  of  the    i         group  to  the  energy 
level  of  1.  44  ev. 

Another  solution  for  such  a  source  distribution  shown  for  the  Pu-Be 
source  in  Fig.    1  would  be  found  from  a  numerical  integration  of  the  point 
kernel  over  the  source  spectrum  (43), 

v    rT  -b.2(t-t    ) 

q(r,  t)  =  2    J0  dTs   Si(Ts)>T"i  <r  )  e  S    .  (n; 

i  =  l 

S.(t)    is  the  amplitude  of  the    i         harmonic  from  a  source  of 
age     t    .       *!'.(  r  )     represents  the  wave  function  satisfying  the  age 
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Fig.  I.   Calculated  and  measured  neutron  energy  spectrum  for 
a    Pu_Be  neutron  source. 
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equation  and  the  boundary  conditions,     and    B.2    represents  the  buckling 
of  the  distribution    $.(  r  ). 

An  approximate  fit  of  the  quantity  proportional  to  the  slowing  down 
density  at  indium  resonance  energy  may  be  made  with  a  finite  number 
of  the  terms  in  Eq.    12. 

f.   C.e  (14) 

i  =  l 

In  terms  of  the  experimental  distribution,   Eq.    14  is 

A  =  V  B.  e  '  (15) 

res      Lj      i 

i  =  l 

The    f.     may  be  found  from  the  experimentally  determined 

amplitudes.     B. ,     and  the  constants    C.     may  be  determined  by 
*  *  '        Ny  y 

Eq.    10.     B.     is  proportional  to    f.  C.     and     /       f .    =  1 .     The  combi- 

l  11  Lj        i 

nation  of  these  relationships  gives 


f;     ■     -N1-2 W 


B  ,/C. 
r 

N 

The  result  of  this  analysis  gives  an  analytical  expression, 
Eq.    14,    for  the  absolute  slowing  down  density  to  indium  resonance 
energy  if  the  absolute  source  strength    Q0    is  known. 

An  experimentally  determined  Gaussian  range  to  indium  reson- 
ance energy  may  be  extended  to  represent  the  Gaussian  range  to  thermal 
energy  by  use  of  the  definition  of  the  age,    Eq.    5,    and  by  use  of  the  re- 
lationship between  age  and  Gaussian  range,    Eq.    9. 
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T<E01-V  =  T(Eoi'EIn>  +  T(EIn'    Eth' 


r2  -2 

rth-i 


6  2, 


D  d_E_ 

E 


(17) 


(18) 


Using  the  Gaussian  ranges  to  thermal  energy,   the  absolute 
expression  for  slowing  down  density,    Eq.    14,    gives  the  number 
of  neutrons  slowing  down  to  thermal  energy  per  unit  volume  per 
second,   i.e., 

7r 


^k-i 


i  =  l 


This  eM  i.ipolation  of  the   resonance  energy  Gaussian  ranges 
to  thermal  energy  assumes  that  the  moderating  process  described 
by  the  agi:  equation  applies  until  the  neutron  energy  corresponds  to 
thermal  energy.     At  this  point  slowing  down  stops  and  the  neutron 
is  considered  to  diffuse  at  a  constant  thermal  energy. 

The  solution  to  the  thermal  diffusion  equation  in  the  infinite 
column  of  rectangular  sides    a    and    b,   using  as  a  source  term 
Eq.    19,    the  expression  for  absolute  slowing  down  density  to  thermal 
energy,   is  derived  in  Appendix  A  and  given  here  as 
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where 


and 
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It  is  evident  that  if  the  nuclear  constants  of  the  moderator  are 
known,    the  absolute  neutron  flux  along  the     Z   axis  may  be  defined 
from  Eq.    20.     Once  this  absolute  flux  is  defined,   various  thermal 
neutron  detectors  can  be  standardized  against  it  so  that  a  normali- 
zation factor  will  convert  from  count  rate  to  flux.     Such  detectors 
will  be  referred  to  as  standardized  detectors. 

To  correct  for  the  fact  that  neutron  measurements  cannot 
actually  be  made  in  a  column  of  infinite  height,   the  experimental 
thermal  measurements  must  be  adjusted  so  that  it  appears  that  they 
were  made  in  a  pile  of  infinite  height.     Actually  the  neutron  flux  goes 
to  zero  at  the  upper  extrapolated  boundary,     h,     of  the  actual  finite 
column.     Actual  pile  measurements  may  be  adjusted  by  the  familiar 
fundamental  mode  end  correction  factor  used  in  diffusion  length  measure- 
ments (11), 


1 


-2(h-z)/b 


I 


(23) 


The  solution  to  the  diffusion  equation  in  a  finite  medium  of  sides 
2a,    2b,    and  2c,   using  as  a  source  term  the  distributed  thermal  source 
given  by  Eq.    6,    is  derived  in  Appendix  B  and  given  here  as 
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In  this  solution    Fm(x),    Fm(x0),    Gn(y),    Gn(yo).    Hq(z),    Hq(zq^ 

are  given  in  Eqs.    7. 

If  the  thermal  diffusion  equation  in  this  geometry  is  sourced  by 
a  thermal  point  source  the  solution,   derived  in  Appendix  C,    is 
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In  this  equation,  \.  represents  the  first  collision  mean  free  path 
of  the  neutrons  in  the  i  energy  group.  \.  may  be  determined 
according  to  an  equation  given  on  page  200  of  Reference  43. 


(25) 


Some  note  should  be  made  of  the  inaccuracy  in  the  assumption 

of  a  physical  point  neutron  source.      The  source  to  the  age  equation 

should  actually  be  given  by  the  distribution  of  first  collisions  (27). 

It  would  be  introduced  into  the  age  equation  as  a  function  of  the  mean 

free  path  and  the  distance  from  the  source  location 
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EXPERIMENTAL   FACILITIES 
General  Pile  Description 

The  Kansas  State  University  Standard  Graphite  Pile,    shown  in 
Fig.    2,   in  which  all  measurements  were  made,   consisted  of  a  rectang- 
ular parallelepiped,     68  inches  square  and  100  inches  high,    resting  on 
a  concrete  foundation.      The  column  was  constructed  by  stacking  layers 
of  machined  reactor  grade  graphite  blocks,    approximately  4  inches  in 
cross  section  and  of  various  lengths.     In  stacking,   the  long  dimension 
of  the  graphite  blocks  was  altered  by  90°  cross-lapping  from  layer  to 
layer.     As  seen  in  Fig.    2,    certain  of  the  graphite  blocks  were  drilled 
to  a  diameter  of  1.75  inches  along  the  lengths  of  the  blocks.     For  this 
work,   in  the  effort  to  make  a  homogeneous,    solid  moderator,    each  of 
the  fuel  ports  contained  graphite  cylinders  measuring  1.625  inches  in 
diameter  with  each  being  22.68  inches  in  length,   and  of  the  same 
material  as  the  graphite  blocks. 

The  graphite  cylinders  rested  on  the  bottoms  of  the  fuel  ports,    leav- 
ing crescent  shaped  air  gaps  between  the  tops  of  the  cylinders  and  the 
tops  of  the  fuel  ports.     The  distance  between  the  tops  of  the  cylinders 
and  the  tops  of  the  fuel  ports  was  0.  125  inches.      The  blocks  along  the 
central  vertical  axis  of  the  column  contained  horizontal  foil  slots  with 
a  cross  section  of  1.281  inches  by  0.343  inches.      These  slots  were  such 
that  graphite  foil  stringers,    shown  in  Fig.    3,    containing  indium  foils 
with  their  covering  material,    could  be  placed  at  approximately  4  -inch 
intervals  up  the  central  axis  of  the  column.     The  density  of  the  solid 


Fig. 2.    KSU   Standard    Graphite  Pile. 
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graphite  blocks  was   1.683  gm/cm3.     Steichen  (39)  calculated  the  per- 
centage of  air  voids  to  be  approximately  0.  3  per  cent.     The  resulting  ef- 
fective  density  of  the  pile  was  1.678  gm/cm3. 

A  summary  of  the  physical  specifications  of  the  one  curie  Pu-Be 
neutron  sources  is  given  in  Table   1. 

A  source  was  held  in  a  graphite  positioning  cylinder,    shown  in 
Fig.    3,    of  inner  diameter  such  that  the  sources  fit  snugly  in  the  cylin- 
ders and  of  outer  diameter  such  that  the  cylinders  fit  closely  in  the 
fuel  ports  of  the  graphite  pile. 

Standard  Pile  Configuration 

For  the  standard  pile  analysis,    a  source  in  its  graphite  positioning 
cylinder,   was  located  in  the  column  in  the  bottom  row  of  fuel  ports 
(just  above  the  pedestal)  such  that  the  central  axis  along  the  length  of  the 
source  containers  coincided  with  the  horizontal  central  axis  of  the  fuel 
port,   and  also  so  that  the  center  of  the  source  coincided  with  the  verti- 
cal central  axis  of  the  pile.     In  this  position,  the  center  point  of  the 
source  arrangement  was  26  inches  from  the  bottom,    74  inches  from  the 
top,   and  34  inches  from  any  side  of  the  column. 

A  schematic  diagram  of  the  standard  pile  showing  the  various 
detector  locations  in  relation  to  the  source  is  shown  in  Fig.    4. 

Finite  Pile  Configuration 

For  the  finite  pile  analysis,    a  source  was  positioned  in  the  plane  of 
the  fifth  layer  of  fuel  ports,   in  the  same  centering  arrangement  within  the 
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Table  1.     Summary  of  physical  specifications  of  Plu- 
tonium-Beryllium neutron  sources  (8). 


Source 
number 

Grams  Pu 

Grams  Be 

Neutron  emission 
rate  (n/sec) 

365 

7.87 

15.99 

1.54  x  106 

366 

7.87 

16.01 

1.61  x  106 

367 

7.86 

15.89 

1.71  x  106 

368 

7.86 

15.88 

1.58  x  106 

369 

7.86 

16.09 

1.60  x  106 

The  inner  container  of  each  source  is  0.85  inches  in  diameter, 
0.90  inches  in  height,    and  is  made  of  tantalum.     The  outside  con- 
tainer of  each  source  is   1 .  02  inches  in  diameter ,    1.30  inches  in 
height,    and  is  made  of  stainless  steel.     All  sources  are  sealed  by 
welding.     The  melting  point  of  tantalum  is  5,425°  F.     The  melting 
point  of  18-8  stainless  steel  i.s   2,600°  F. 

The  sources  were  calibrated  by  comparison  to  within  ±2  per 
cent  of  standards  calibrated  at  the  Los  Alamos  Scientific  Laboratory. 
The  absolute  accuracy  of  the  Los  Alamos  standard  is  reported  as 
±5  per  cent  (7),     thus  giving  an  absolute  accuracy  of  ±  7  per  cent  to 
the  sources  listed  above. 
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Fig. 4.     Schematic    diagram   of    KSU  Standard  Graphite   Pile. 
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plane  as  was  used  with  the  standard  pile  arrangement.     Cadmium  sheets, 
approximately  0.  020  inches  thick  were  inserted  in  the  plane  separating 
the  pedestal  from  the  first  row  of  fuel  port  blocks  in  order  to  create  an 
artificial  boundary  at  that  point.     Provision  for  this  arrangement  was 
made  at  the  time  of  the  original  construction  of  the  pile  by  inserting 
strips  of  1 1  00  aluminum  alloy  having  a  width  of  0.  50  inch  and  a  thick- 
ness of  0.032  inch.     The  aluminum  strips  were  spaced  such  that  sheets 
of  cadmium  approximately  8  inches  wide  could  be  inserted  easily. 

The  two  uppermost  layers  of  graphite  were  removed  from  the  pile 
for  the  finite  geometry  measurements.     This  left  the  pile  in  such  a 
state  that  it  approximated  a  cube  with  sides  of  68  inches,    and  with  one 
source  in  the  center  of  this  cube.     A  schematic  view  of  this  pile  arrange- 
ment is  shown  in    Fig.    5. 

Foil  Counting  Facilities 

The  indium  foils,    shown  in  Fig.    6,   were  1.00  inches  in  diameter, 
had  a  thickness  of  approximately  0.  005  inches,    and  weighed  an  average 
of  0.4592  gm,    giving  a  mg/cm2  thickness  of  92.7.     The  foils  were  chosen 
such  that  their  weights  did  not  differ  by  more  than  one  per  cent. 

The  cadmium  covers,    shown  in  Fig.    6,   were  constructed  by  press- 
ing them  from  a  cadmium  sheet.     The  covers  had  an  inside  diameter  of 
slightly  more  than  one  inch  and  a  wall  thickness  ranging  from  0.  021 
inches  to  0.022  inches.     In  addition  to  the  cover,   the  foils  were  covered 
on  both  sides  with  cadmium  disks  of  the  same  thickness  as  the  covers, 
giving  an  overall  average  cadmium  covering  thickness  of  0.043  inches. 
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Fig.  5.    Schematic  diagram   of   finite   pile. 
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Aluminum  covers  were  made  in  the  same  manner  for  the  purpose 
of  protecting  and  supporting  the  foil  during  bare  foil  measurements. 
Aluminum  disks  were  inserted  between  the  foil  and  the  cover  in  order 
to  duplicate  the  irradiation  geometry  used  in  the  cadmium  covered 
foil  measurements. 

The  foil  stringers,    0.125  inches  by  0.325  inches  by  12  inches, 
shown  in  Fig.    3,   were  constructed  from  the  same  material  as  the 
graphite  blocks.      The  stringers  were  recessed  so  that  the  foils  with 
their  covering  material  could  be  positively  positioned  in  the  center  of 
each  stringer. 

The  foil  counting  system  shown  in  Fig.    7  consisted  of  a  B.    J. 
Electronics     Model  DD7     continuous  flow  proportional  counter,    a  Baird 
Atomic     Model  255     proportional  counter  preamplifier,    a  Baird  Atomic 
Model  132     scaler,   and  a  Baird  Atomic     Model  323     timer.     The  KSU 
Nuclear  Engineering  Inventory  No.    of  the  scaler  was   145,    and  the  pro- 
portional counter,   369.     The  local  identification  of  the  preamplifier 
was  KSU  Nuclear  Laboratory  Inventory  No.    10025.     The  proportional 
counter  had  a  tungsten  collector  wire  with  a  one  half  inch  loop.     The 
counting  gas  was  Olin-Matheson  P-10,    a  mixture  of  10  per  cent  methane 
and  90  per  cent  argon. 

The  proportional  counter  was  housed  in  a  lead  brick  shield  such  that 
the  minimum  shield  thickness  at  any  point  was  at  least  2  inches.     A  0.  020 
inch  thick  aluminum  tray  was  constructed  to  facilitate  the  removal  and  the 
positioning  of  the  foil  in  the  counter  well. 

The  major  portion  of  the  foil  counting  system  just  described  was  due 
to  the  efforts  of  Steichen  (39). 
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Neutron  Probe  Facilities 

The  counting  system  used  in  the  neutron  probe  measurements  is 
shown  in  Fig.    8.     Measurements  were  made  with  both  a  BF3  probe, 
Nuclear  Engineering  Inventory  No.    97,    and  a  Scintillation  probe, 
N.   E.   No.    99.     The  BF3  probe     counting  system  consisted  of  a  Nuclear 
Chicago     Model  NC  - 202     BFj  neutron  probe,    a  Nuclear  Chicago 
Model  1062     preamplifier,    a  Baird  Atomic     Model  132     scaler,    a 
John  Fluke     Model  400  BDA     high  voltage  power  supply,    a  B.    J. 
Electronics     Model  DM1 -D     count  rate  meter,    and  an  Esterline  Angus 
Model  A.W.    strip  chart  recorder. 

The  scintillation  probe  counting  system  was  identical  except  that 
a  Nuclear  Chicago,    Model  DS8-10,    scintillation  detector,   was  sub- 
stituted for  the  NC-202  BF3  probe  and  its  Model  1062  preamplifier. 

The  basic  detector  portion  of  the  DS8-10  scintillation  probe  had 
a  transitorized  preamplifier  and  a   10-stage  photomultiplier  with  cesium 
antimony  dynodes,   housed  in  a  container  9-1/4  inches  long  by  1  inch 
diameter.     The  SN-6  neutron  probe  portion  had  a  thermal  neutron 
crystal  which  was  a  phosphor  of  silver  activated  zinc  sulfide,    95  per 
cent  enriched  with  boron- 10.     It  was  attached  at  the  end  of  an  aluminum 
covered  lucite  light  pipe,    42  inches  long  by  5/16  inch  outside  diameter. 

The  BF3  probe  active  volume,   with  length  of  l/2  inch  and  diameter 
of  3/l6  inch,    contained  B10  F3  gas  at  a  pressure  of  70  cm  Hg  and  at  an 
enrichment  of  96  per  cent.     The  Model  1062  transistorized  preamplifier 
had  a  length  of  6.25  inches  and  a  diameter  of  1 .  5  inches. 
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The  BF3  probe  is  shown  with  the  counting  system  in  Fig.    8;     the 
scintillation  probe  is  shown  in  Fig.    12. 

Traversing  Mechanism 

For  the  purpose  of  supporting  these  two  detectors  in  a  variety  of 
positions  in  the  pile,    a  traversing  mechanism  was  designed  and  built. 
A  drawing  of  this  mechanism  is  displayed  as  Fig.    9.     Figures  10  and 
11  include  photographs  of  the  complete  assembly,    and  Fig.    12  shows 
a  close-up  view  of  the  positioning  device  of  the  mechanism. 

The  mechanism  was  designed  such  that  the  detector  could  be 
positioned  positively  at  any  elevation  between  the  limits  of  19  inches 
and  101  inches.     The  detector  could  be  positioned  horizontally  into 
the  pile  either  by  hand  adjustment  or  motor  powered  adjustment.     A 
cam  action  engaging  device  allowed  the  detector  to  be  moved  horizon- 
tally at  the  rate  of  one  inch  per  minute  by  the  motor  and  lead  screw, 
or  it  allowed  the  detector  to  be  positioned  freely  by  hand.      The  mechan- 
ism was  designed  such  that  a  maximum  horizontal  traverse  of  40  inches 
was  available. 

The  complete  traversing  mechanism  could  be  moved  to  a  different 
horizontal  position  on  the  face  of  the  pile  as  a  result  of  caster  wheels 
mounted  on  the  base.     The  base  of  the  mechanism  could  be  positioned 
firmly  in  a  desired  position  with  respect  to  the  floor  by  use  of  the  level- 
ing screws  at  each  corner  of  the  base,   which  allowed  the  caster  wheels 
to  be  raised  from  the  floor. 


29 


30 


MTTECTOR 


MSCH1PTIQN 

FORMICA       PULLEY     89  126    4  ^   STL.    SHAFT 

PULLEY      SUPPOPT,     Z;r3'4*Vfe     STL.    BAR 

Mb     5.S.     CABLC 

COUNTER      WEIGHTS 

TOP     BEAM,    2*2* '-'ft    ANGLE 

TUBE"     HOLDER, TOP,    2    SCH.     ©O    BLACK    PIPE 

VERTICAL   TUBES,    2.004    O.O.     S.S. 

VERTICAL    BRACE,      I      BLACK      IRON      PIPE 

END     BLOCK,     2  VZ  *  I      AL      BAR 

SLIDING    BLOCK,    2  >  I     AL.    BAR 

GUIDE    TUBES,    1.378    O.D.    S.S. 

CLAMP    rfTIGHTNER,    2      SCH.      BO    BLACK   PIPE 

DETECTOR      HOLDER,     FISHER    CLAMP 

DETECTOR 

LEAD    SCREW,    34-IONC 

LEAD     SCREW     BEARINGS,   BOSTON     NO.  1 6(6 

SPACER,    I  '-i    BLACK    PIPE 

MOTOR,     BOD1NE;     NO.     BSI30E-I60M 

MICROSWITCH,  NC     5A,      I23"VA.C 

>4  A.,  SLQ-BLO,     FUSE     4    FUSE     HOLDER 

MOTOR    SWITCH,    S.RD.T      SA.     I25V.AC 

MOTOR     CAPACITOR,    1.2      MFD,    22QVA.C 

THREADED     5PACER    AND     CAM   CLAMP 

CONTROL     MOUNT,  2x2  » '/a   ANGLE 

6  FT        STL      TAPE 

BUSHINGS,  CLEVE.    GRAPH     NO.    5"3I-C 

GEAR,   BOSTON     NO.    Y4B72 

LEVEL     INDICATOR,     5EARS  *  ROEBUCK 

BASE",  A    PCS.    3x3xW  ANGLE,  S   PCS   '4  PLATE  i^ 

CASTE R,COLr>ON     NO.  5329I-HARD,  5S212-SOFT     w\ 

TUBE     HOLDER, BOTTOM,    Z    BLACK    PIPE 

ANCHOR,    2»'/4      IRON      STRAP 

LEVELING      SCREW,    5X*     MACHINE    &OLT 


RIGHT    SIDE 


0 


<§K 


®- 


w 


@^ 


@, 


FRONT 


RIGHT    SIDE 


Fig.  9.  Traversing  mechanism  for  graphite  pile,  assembly  drawing. 


Fig.  10.  Traversing   mechanism  and   standard   pile. 


Fig.  II.   Traversing   mechanism    and    finite   pile. 
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In  order  to  accommodate  these  neutron  probes  in  a  maximuni 
number  of  positions  and  with  the  minimum  amount  of  created  void 
space,  the  three  1.652  inch  diameter  graphite  cylinders,    shown  in 
Fig.    13,  were  drilled  along  their  central  axis  to  a  diameter  of  3/8-inch. 
Their  use  to  accommodate  a  probe  may  be  seen  in  Fig.    11. 

EXPERIMENTAL   PROCEDURE 
Foil  Measurements 

All  foil  experiments  were  performed  using  one  of  the  neutron  sources 
described  in  Table  1.     The  source  was  placed  in  position  in  the  pile  as 
described  in  the  section  on  experimental  facilities.     The  foils  in  their 
appropriate  covering  material  were  placed  in  the  foil  stringer.     The 
foil  stringers  were  then  placed  in  the  pile  such  that  the  indium  foils  were 
centered  on  a  line  along  the  vertical  axis  of  the  column,   directly  above 
the  center  of  the  neutron  source. 

The  foils  were  allowed  to  remain  in  the  pile  at  least  7  hours  until 
they  reached  equilibrium  saturation  activity  with  respect  to  the  54  min- 
ute   half-life  isotope.     The  foil  nearest  to  the  source  was  removed  and 
a  spare  foil  in  an  identical  covering  material  was  placed  in  the  vacated 
position.     The  removed  foil  was  allowed  to    "cool"  for  3  minutes  to 
assure  the  decay  of  the  13  second  half-life  isotope,    and  was  then  counted 
for  a  period  of  27  minutes.     At  the  end  of  the  counting  period,   the  foil 
second  closest  to  the  source  was  removed  and  the  foil  which  had  just 
been  counted  was  immediately  placed  in  the  vacated  position.     This 
procedure  was  followed  until  all  foils  in  the  pile  had  been  counted.     A 


Fig.  13.    Special   graphite     cylinders   for  accommodation  of 
neutron    probes. 


total  of  eleven  sets  of  data  using  cadmium  covered  indium  foils  were 
taken  using  this  procedure.     Of  these  eleven  trials,    5  were  made  count- 
ing the  activity  of  the  source  side  of  the  foil,    and  6  were  made  count- 
ing the  side  of  the  foil  away  from  the  source.     A  total  of  4  sets  of  data 
were  taken  using  aluminum  covered  indium  foils.     Of  these  4  trials. 
2  were  made  on  each  side  of  the  foil.      The  background  count  was  taken 
for  a  period  of  27  minutes  before  and  after  the  loimliny  periods  for 
each  of  the  sets  of  data  .     The  dala  (tu-  cadmium  covered  indium  foils 
were  taken  by  Steichen  (39). 

The  counter  was  operated  at  1  ,  950  volts,    and  with  a  pulse  height 
sensitivity  setting  of  0.  2.     A  number  of  trials  were  made    by  Steichen 
using  various  sensitivity  settings,    foil-collector  wire  spacings,    and 
sizes  of  collector  wire  loops,    to  obtain  an  optimum  operating  plateau 
and  minimum  counter  background. 

The  counting  system  was  checked  for  stability  before  and  after 
each  trial  with  standard  Ra-D-E  source,    Nuclear-Chicago  Serial 
No.    2354.     All  foil  data  used  in  pile  standardization  were  taken  with 
the  counting  system  in  such  a  condition  that  the  count  rate  from  the 
standard  Ra-D-E  source  was  within  200  counts  of  17,  048  counts  per 
minute. 

Neutron  Probe  Measurements  in  Standard  Pile 

All  experiments  with  the  neutron  probes  were  performed  using  one 
of  the  neutron  sources  described  in  Table   1.     Probe  measurements  were 
made  in  both  the  foil  slots  and  the  fuel  ports.     Measurements  in  the  fuel 
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ports  were  made  using  the  special  graphite  cylinders  described  in  the 
section  entitled  "EXPERIMENTAL   FACILITIES.  " 

Probe  data  for  the  standard  pile  analysis  were  taken  such  that 
10,000  counts  were  obtained  in  each  selected  position.     For  the  BF3 
probe,   this  was  done  with  7  series  of  runs  getting  1,  000  counts  in  each 
position,    and  one  run  getting  3,  000  counts  in  each  position.     For  the 
scintillation  probe,    this  was  done  with  a  series  of  10  runs,    getting 
1,  000  in  each  selected  position. 

BF3  probe  data  were  taken  at  an  operating  voltage  of  1,450  volts, 
and  a  pulse  height  sensitivity  setting  of  0.8.     Scintillation  probe  data 
was  taken  at  the  same  pulse  height  sensitivity  but  with  an  operating 
voltage  of  11,  000  volts. 

Before  each  series  of  runs,    10  one-minute  counts  were  taken  in  a 
standard  position  to  check  for  proper  distribution  and  counter  repro- 
ducibility.    All  BF3  probe  measurements  were  taken  at  an  operating 
condition  such  that  a  count  rate  of  1,  764    ±  25  counts  per  minute  was 
obtained  at  a  position  15.82  cm  directly  above  source  number  368  in 
the  standard  pile  configuration.     All  scintillation  probe  measurements 
were  taken  at  an  operating  condition  such  that  a  count  rate  of  2533  ±57 
counts  per  minute  was  obtained  in  this  same  standard  position.'     The 
variance  reported  on  these  count  ranges  is  the  standard  deviation  of  a 
single  observation. 

The  count  rate  meter  and  the  strip  chart  recorder  were  used  to 
monitor  the  counting  during  periods  of  operator  absence,    thus  insuring 
that  no  erratic  data  were  recorded. 
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Neutron  Probe  Measurements  in  Finite  PiK- 

Al!  measurements  for  the  experiment  in  the  finite  geometry  which 
npared  the  various  source  conditions  to  the  thermal  diffusion  equa- 
tion were  made  with  the  standardized    BF3    probe  operaimg  at  the  con- 
ditions described  in  the  previous  section. 

The  medium  was  mapped  by  means  of  the  special  graphite  cylinders 
placed  in  the  five  hole  by  five  hole  array  of  fuel  ports  of  the  upper, 
northwest  l/8th  of  the  finite  medium  (see  Fig.    5).     By  selecting  posi- 
tions at  4-inch  intervals  along  the  horizontal  traverse  into  the  special 
graphite  cylinders,    9  positions  were  available  in  each  of  24  fuel  ports, 
and  8  positions  were  available  in  the  one  fuel  port  containing  the  source, 
thus  making  a  total  of  224  positions  to  be  counted  for  each  pile  mapping. 
This  procedure  resulted  in  a  3-dimensional  mapping  using  a  net  dimen- 
sion of  4  inches  by  8  inches  by  8  inches. 

The  data  were  taken  such  that  in  each  position,    either  1,  000  counts 
were  obtained  or  10  minutes  counting  time  had  elapsed,   depending  on 
which  condition  was  met  first.     Each  selected  position  was  counted  for 
at  least  2  minutes.     Four  complete  traverses  of  the  224  positions  were 
made  to  get  the  average  at  each  position  as  reported  in  the  results. 

DATA   PRESENTATION   AND   ANALYSIS 

Pile  Standardization  Data  Treatment 

The  raw  data  averages  used  in  the  standard  pile  analysis  were  cor- 
rected for  background  and  tabulated  in  Table  2.     All  data  used  in  the  pile 


39 


standardization  were  based  on  source  No.    368  used  by  Steichen,   having 

an  emission  rate  of  1.58  x  1  0    neutrons  per  second  (see  Table   1).     It 

is  important  to  note  that    A_    ,    the  activity  of  an  aluminum  covered  or 

bare  indium  foil  and    A_  ]T  ,     the  activity  of  a  cadmium  covered 

CdlnCd  ' 

indium  foil,   has  always  been  reported  in  this  work  as  the  total  counts 
over  the  27  minute  counting  period. 

The  foil  data  may  be  adjusted  to  give  activity  at  the  time  of  removal 
by  well  known  decay  laws  which  results  in  the  equation, 


Total  counts  in  _„ 

,-       .       .  30  min  . 

27  minute  counting  r  ~\q*- 

period  =      J  Agat    e      d     dt.  (27) 

3  min 

Using  X.  as  0.  012815  min  for  the  appropriate  decay  constant 
of  indium,  the  saturation  activity  at  the  time  of  removal  from  the 
pile  may  be  found  from 

=    (    T°tal  c°u"ts  in  Z1  ™in"   )    (  0.0422)   .  (28) 

Sat  ute  counting  period  '    x  '  ' 

This  may  be  considered  to  be  the  saturation  activity  since  the  foils 
were  always  irradiated  for  seven    half-lives  or  more. 

Column  2  of  Table  2  lists  the  cadmium  covered  indium  foil 
data  as  measured  by  Steichen  in  his  determination  of  the  age  of 
neutrons  from  a  Pu-Be   source  in  graphite  (39).     Columns  4  and  5 
list  the  counts  per  minute  obtained  from    BF3    probe  measurements 
and  scintillation  probe  measurements.     The  error  quoted  in  all  cases 
is  the  standard  deviation  of  the  mean.     It  includes  only  the  statistical 


Table  2.     Summary  of  standard  raw  data. 


40 


z 

ACdInCd 

A 
In 

ABF3 

A_ 

Sc 

cm 

total  counts 
in  27  minute 
ount  period 

total  counts 
in  27  minute 
count  period 

counts  per 
minute 

counts  per 
minute 

5.63 

35,199.0  ± 

92.0 

201,715  ±3,367 

1,856.00  ±14.0 

2,823.0  ±19.0 

15.82 

27,668.0  ±101.0 

184,099  ±     604 

1,764.00  ±   8.0 

2,501.0  ±18.0 

25.96 

18,543.0  ± 

42.0 

152,790  ±    368 

1,528.00  ±14.0 

2,217.0  ±25.0 

36.11 

10,752.0  ± 

34.0 

1,228.00  ±10.0 

1,716.0  ±14.0 

46.31 

5,686.0  ± 

35.0 

89,840  ±    317 

973.00  ±   8.2 

1,341.0 ±  16.0 

56.45 

2,558.0  ± 

12.0 

62,977  ±1,057 

700.00  ±  6.1 

948.9  ±  8.7 

66.15 

1,174.0  ± 

3.0 

515.90  ±  6.9 

692.0  ±6.1 

76.78 

468.5  * 

6.4 

32,437  ±     286 

371.60  ±2.5 

509.6  ±4.2 

86.94 

209.8  ± 

6.4 

23,661  ±     116 

272.90  ±  2.5 

374.2  ±4.3 

97.12 

93.1  ± 

6.1 

196.40  ±2.2 

256.7  ±  2.5 

117.U 

8,763  ±     148 

100.40  ±  2.0 

132.4  ±1.8 

137.76 

4,287  ±      95 

51.30 

66.5  ±  1.1 

158.08 

24.63 

32.0  ±  0.9 

178.40 

7.54 
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error  of  the  experimental  data.     Steichen's  data  for    A_  ,.    _  ,    were 
a  combination  of  data  taken  with  four  sources  and  data  taken  with  one 
source.     The  corrected  activities  for  four  source  geometry  were 
normalized  to  the  one  source  geometry  at  46.31  cm.     Steichen  felt 
that  the  four  source  data  gave  better  statistics  than  the  one  source 
data  at  points  far  from  the  source.     At  these  points,    the  correction 
for  the  finite  size  of  the  four  sources  became  small. 

The  purpose  of  the  foil  measurements  was  to  get  quantities  pro- 
portional to  the  thermal  flux  and  the  slowing  down  density  at  indium 
resonance  energy.     The  activity  due  to  thermal  neutrons  only  was 
determined  by  means  of  a  cadmium  difference  method.     The  foil  was 
exposed  to  neutrons  at  a  given  position,    first  uncovered  and  then 
covered  with  a  layer  of  cadmium  0.  043  inches  thick.     Since  cadmium 
absorbs  the  thermal  neutrons,   the  activation  of  the  cadmium  covered 
indium  foil  was  due  primarily  to  epithermal  neutrons  while  the  un- 
covered indium  foil  activation  was  due  to  both  thermal  and  epithermal 
neutrons.     The  difference  in  these  two  activations,   the  cadmium  dif- 
ference,   after  account  was  taken  of  the  epithermal  neutron  absorption 
in  the  cadmium  cover  by  the  factor    F    ,     (40),   was  taken  to  be  pro- 
portional to  the  thermal  flux. 

<"V>th  =    Kth    Ath     "    K°th    F'  (Aln  -  FCdACdInCd  )  <29> 

In  this  expression,     K        is  the  constant  of  proportionality 
between  flux  and  cadmium  difference  activation  which  standardizes 
the  detector.       F1    is  a  term  containing  considerations  of  perturba- 
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tions  and  corrections  introduced  through  experimental  limitations  and 
finite  foils.     K°      is  the  constant  of  proportionality  which  is  independent 
of  these  corrections  and  perturbations. 

The  term    F1    may  include  the  following  considerations: 

a)  activation  of  indium  by  higher  energy 

resonance  neutrons 

b)  absorption  of  epithermal  neutrons  in  the 

cadmium  filter 

c)  finite  size  of  source  and  detector 

d)  difference  in  foil  weights  and  size 

e)  unequal  activation  of  two  sides  of  foil 

f)  flux  depression  in  vicinity  of  foil 

g)  self-shielding  of  the  center  of  the  foil 

by  its  outer  layers 

h)    mutual  shadowing  of  neighboring  foils 

i)    hardening  of  the  neutron  spectrum 
through  the  foil 

j)    obliquity  of  foil  with  respect  to  the  source. 

These  correction  factors  are  discussed  further  in  Appendix  H. 
Reasons  for  consideration  or  rejection  of  individual  correction  factors 
are  given,   and  the  correction  factors  (a)  through  (e)  which  were  used 
in  this  work  are  discussed  quantitatively.     The  standardization  con- 
stants determined  in  this  thesis  may  be  used  only  in  graphite  since 
many  of  the  neglected  correction  factors  depend  upon  the  diffusing  medium. 

As  outlined  in  the  theory,   the  slowing  down  density  at  indium  reson- 
ance energy  may  be  taken  to  be  proportional  to  the  activation  of  a  cadmium 
covered  indium  foil. 
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q  =  K  A  =  K°         F'  A_  ,T    _,  (30) 

^res  res      res  res  CdlnCd  '      ' 

F1    again  represents  appropriate  correction  factors  that  may  be 
applied  to  the  foil  activations. 

It  must  be  kept  in  mind  that  for  a     l/v    absorber  such  as  indium 
in  the  thermal  range  and  for  a  resonance  absorber  in  the  epithermal 
range,   the  foil  activation  is  more  correctly  proportional  to  the 
density  of  thermal  neutrons  or  to  the  density  of  resonance  neutrons 
(16).     Hence,    for  a  foil  of    l/v    cross  section,     A        is  proportional 
to    n     ,     not    (nv)  ,  ,     making    K        depend  upon  the  energy  distribu- 
tion of  the  thermal  neutrons.     Thus,   foils  standardized  in  graphite 
for  thermal  flux  cannot  readily  be  extrapolated  for  use  in  any  other 
medium  where  the  energy  distribution  of  the  thermal  neutrons  may 
be  different. 

In  the  case  of  a  slowing  down  density  standardization,    where 

A^»c    is  proportional  to    n  and    q         ,     the  standardization  con- 

res  res  res 

stant  for  converting    A  to    q  cannot  be  directly  used  in 

°        res  ^res  ' 

media  other  than  graphite  (16)  since 

Ires  "     ^r"3     ^     ~-    Kres  Ares  <31> 

where    \    represents  the  scattering  mean  free  path.     To  convert 
from  graphite  to  water,   the  slowing  down  density  standardization 
constant  in  water  would  have  to  be  expressed  (16)  as 

K  (water)  =  K  (graphite)    X  (graPhite)£  (water)  ,      . 

res  res  iB      r         '    ,    ,  .   „  .  ,  '■I 

\  (water)  §  (graphite) 
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Table   3  presents  the  corrected  foil  activations     A  and    A  ^ 

which  are  considered  to  be  proportional  to  the  slowing  down  density 
at  indium  resonance  energy  and  the  thermal  neutron  flux,    respectively. 
The  end  correction  factor  given  by  Eq.    23  in  the  theory,   which  adjusts 
the  finite  pile  measurements  to  make  them  appear  as  if  they  had  been 
taken  in  an  infinite  pile,    is  given  in  Column  4  of  Table  3.     Columns  5, 
6,    and  7  give  the  adjusted    infinite  pile  measurements  with  foils,   the 
BF3    probe  and  the  scintillation  probe.     Perturbations  of  the  thermal 
flux  by  the  two  probes  have  not  been  considered.      The  count  rates 
from  these  two  instruments  have  been  assumed  to  be  proportional  to 
the  thermal  neutron  density,    so  that  they  are  subject  to  the  restriction 
that  the  standardization  constant  found  in  graphite  is  not  useable  In  a 
medium  where  the  energy  distribution  of  the  thermal  neutrons  may  be 
different.     The  thermal  standardization  constants  for  these  detectors 
will  be  found  from  the  relationships 

(™>th  =  Kth<BF>>RBF,  (li) 

and 

("*)*  =  Kth  <Sc>  RSc  (i4) 

where     K     signifies  count    rate. 

Empirical  Fit  of  Indium  Resonance  Datu 

il  may  be  recalled  from  the  theoretical  development  ol"  the 
standard  pile  analysis  that  an  empirical  fit  of  the  quantity  pro- 
portional  to  the  glowing  down  density  at  indium   resonance  was 
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required  in  order  to  express  the  slowing  down  density  In  the  analytical 
form  of  Eq.    14.     An  IBM-650  code  was  developed  for  this  purpose  of 
fitting  the  data  with  a  function  of  Gaussian  form.     It  is  given  in 
Appendix  G. 

This  fit  was  made  by  a  trial  and  error  procedure  which  strived  to 
minimize  the  sum  of  the  weighted  squares  of  the  residuals  between 
the  experimental  values  and  the  values  calculated  from  the  empirical 
expression.     Tables  4  and  5  summarize  these  efforts  to  fit  empirical- 
ly   the  resonance  data.     The  effect  of  various  correction  factors  ap- 
plied to  the  resonance  data,    the  effect  of  the  number  of  terms  used  in 
the  finite  summation,    and  the  effect  of  a  weighting  function  were  studied. 

The  best  fit  may  be  considered  on  the  basis  of  two  different  criteria; 
one,    that  the  sum  of  the  weighted  residuals  squared  be  a  minimum;   or 
two,    that  the  residual  between  the  calculated  and  experimental  value 
at  any  point  should  lie  within  the  statistical  standard  deviation  of  the 
experimental  measurement.     A  weighting  function  of  unity  allows  the 
former  criterion  to  be  satisfied  more  satisfactorily.     Scarborough  (33) 
points  out  that  the  residuals  squared  are  best  weighted  with  the  recip- 
rocal of  the  standard  deviation  squared      l/if*     This  weighting  function 
allows  the  latter  criterion  to  be  satisfied  more  completely. 

A  consideration  of  interest  was  the  necessary  number  of  terms 
needed  in  Eq.    14  to   fit  adequately  the  resonance  data.     Since  the  fit 
was  made  to  only  ten  data  points,    there  was  some  limitation  on  the 
number  of  terms  which  could  be  used,    each  term  having  two  parameters. 
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A  fit  of  the  data  with  one  term  was  considered  poor  since  the 
residual  was  much  greater  than  the  standard  deviation  at  all  exper- 
imental points.     A  fit  with  two  terms  was  clearly  an  improvement 
since  the  sum  of  the  weighted  squares  of  the  residuals  decreased 
by  over  a  factor  of  ten.     The  weighted  three-term  fit  of  the'data 
corrected  for  geometry,    higher  resonance  absorption  and  foil  weights, 
was  considered  to  be  the  best  fit.     The  residuals  of  this  fit  were 
within  or  near  the  standard  deviation  of  the  experimental  data  for  all 
points  but  one.     The  fit  at  46.  31  cm  was  consistently  poor  no  matter 
what  manner  of  fit  was  used.     It  may  be  seen  from  either  the  individ- 
ual residuals  or  the  sum  of  the  weighted  residuals  squared  that  a  fit 
with  four  terms  did  not  give  an  improvement  over  a  fit  using  three 
terms.     Inspection  of  the  fourth  term  in  the  fit  reveals  the  effort  of 
the  machine  to  suppress  the  fourth  term  by  giving  it  a  small  amplitude 
and  a  slope  tending  towards  zero,    the   slope  being  equal  to  the  reci- 
procal of  the  square  of  the  Gaussian  range  associated  with  this  term. 

Consideration  of  the  difference  in  foil  weights  had  a  noticeable 

improving  effect  upon  the  empirical  fit.     Little  improvement  was  gained 

from  using  the  correction  factors  for  finite  detector  and  source.     The 

use  of  the  correction  factor  for  absorption  of  neutrons  by  the  higher 

resonances  of  indium  markedly  changed  the  various  amplitudes,    B. , 

the  Gaussians  ranges,      r. ,     and  the  fractional  contributions,     f . . 
i  i 

Using  the  parameters  of  the  appointed  best  fit,    the  expression  for 
the  indium  resonance  data  distribution  was 
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-(     z        )2  -(-z    ^ 

*    *£.  fiA  '  ^   t 


A  =11,112    e      27'3'  +     21,756     e      UM  +  695  e      65-44  (35) 

res  ' 


From  Eq.    14,    the  absolute  slowing  down  density  to  indium  resonance 
energy  was 


q  106  -(— t f  -(_-_  )2  -(-^-)2 

-2" =  1.314  e      "-31  +  2.572  e36'94      +0.824  e      "^  (36) 

Uo 


The  IBM-650  program  required  initial  estimates  of  the  fitting 
parameters  which  were  found  from  a    "stripping-off  "  procedure. 
Shown  in  Fig.    14  are  the  three  "stripped-off"  curves  which  have 
been  superimposed  to  give  the  empirical  fit  of  the  experimental 
points. 

DISCUSSION    AND   RESULTS 

Standardization  Constants 

Table  6  lists  the  values  of  absolute  flux  and  slowing  down  density 
in  the  Kansas  State  University  Standard  Pile,   using  a  source  with  an 
emission  rate  of  1.58  x  1  06  neutrons  per  second.     The  absolute  slow- 
ing density  was  determined  from  Eq.    36.     The  absolute  thermal  flux 
was  obtained  from  the  computer  code  described  in  Appendix  D  which 
numerically  evaluated  the  absolute  flux  defined  by  Eq.    20.      This  code 

used  the  empirically  determined  parameters    f.     and    r.     after    r.     had 

ill 

been  extended  to  thermal  energy  by  Eq.    18.     Figures   15  and  16  show 
the  absolute  flux  and  the  absolute  slowing  down  density  to  indium 
resonance  energy.     Superimposed  upon  these  calculated  curves  are 
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Fig.  14.  Indium   resonance  distribution  as  empirically  fit   by  three 
terms    of    Gaussian  form. 
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Fig.  16.     Slowing  down  density   to  indium    resonance    energy 
in    KSU   Standard    Pile. 
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points  showing  absolute  quantities  as  determined  by  detectors  which 
have  been  standardized  according  to  this  thesis. 

Columns  4  through  7  of  Table  6  list  the  standardization  constants 
of  the  various  detectors  as  they  were  determined  in  the  standard  pile 
analysis.     The  average  standardization  constant  is  listed  at  the  bottom 
of  each  column. 

Figures  17  show  the  variation  of  the  measured  standardization 
constants  for  each  detector  as  a  function  of  the  distance  from  the  source. 
Since  the  spectrum  is  harder  near  the  source,    a  non-random  distri- 
bution of  these  constants  would  indicate  a  possible  energy  dependence 
in  the  assumed  thermal  neutron  detectors.     It  can  be  seen  from  Figs.  17 
that  the  distributions  appear  to  be  random  with  respect  to  the  average 

for  all  detectors  except  the  scintillation  probe.     The  standardization  con- 
stant distribution  of  the  scintillation  probe  suggests  that  the  scintillation 
probe  has  an  appreciable  sensitivity  for  neutrons  of  epithermal  energies. 
In  order  to  determine  a  more  reliable  standardization  constant  for  the 
scintillation  probe,    it  would  be  necessary  to  use  a  cadmium  difference 
technique  such  as  that  used  with  the  indium  foils. 

Subsequent  use  of  these  standardization  constants  has  been  made 
with  the  average  standardization  constant. 

The  absolute  flux  or  slowing  down  density  in  a  graphite  medium  may 
be  found  from  a  particular  standardized  detector  by  merely  multiplying 
its  standardization  constant  by  a  properly  measured  quantity.     In  the 
case  of  the  Kansas  State  University    BF3    probe  or  scintillation  probe, 
this  quantity  is  counts  per  minute,    corrected  for  background.     In  the 
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case  of  KSU  standard  iridium  foils,    the  quantity  is  expressed  in  terms 
of  total  counts  over  the  27  minute  counting  period  as  explained  in  the 
experimental  outline. 

From  Eq.    28,   it  was  shown  that  for  foil  measurements  this  pro- 
perly determined  quantity  could  be  the  saturation  activity  at  the  time 
of  foil  removal  from  the  pile  if  the  standardization  constant  is  used  as 
K(foil)/0.0422. 

Before  a  standardization  constant  may  be  used,    it  must  be  verified 
that  the  counting  system  reproducibility  has  not  changed.     This  can  be 
done  by  reproducing  the  counter  stability  check  points  described  in  the 
discussion  of  the  experimental  methods. 

Confidence  Limits  on  Flux  Assignment 

The  deviation  reported  on  each  standardization  constant  is  the 
standard  deviation  of  the  mean  standardization  constant.     This  deviation 
is  a  measure  of  the  accuracy  to  which  the  detector  was   standardized.     It 
is  not  a  measure  of  the  accuracy  of  the  flux  assignment,    since  as  previous- 
ly pointed  out,    the  accuracy  of  the  flux  assignment  depends  upon  the  ac- 
curacy in  the  knowledge  of  the  source  strength  and  particular  diffusion  para- 
meters. 

Neutron  Sources.     The  neutron  sources,   listed  in  Table  1,    were  pre- 
pared and  calibrated  by  the  Mound  Laboratory.     Their  catalog  (8)  states 
that  the  sources  were  calibrated  by  comparison  with  standards  calibrated 
at  the  Los  Alamos  Scientific  Laboratory.     Since  the  Los  Alamos  source 
has  a  relative  absolute  calibration  ratio  of  1.016,    (7),    compared  with 
the  National  Bureau  of  Standards  source,     a  KSU  source  strength 
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should  be  divided  by  1.  016.     Table   1   explains  the  assignment  of  the  con- 
fidence limit  of   ±  7  per  cent  to  KSU  sources. 

The  Ra-Be  (or,  n)  source  number  38  at  Argonne  National  Laboratories 
has  a  relative  absolute  calibration  ratio  of  1.049  compared  with  the 
Bureau  of  Standards   source  (7).      To  compare  the  KSU  source  with  the 
Argonne  source,   the  KSU  source  should  be  divided  by    1.016/1.049. 
The  Hanford  source  is  not  compared  with  NBS-I  in  Reference  (7). 

Diffusion  Coefficient.       A  diffusion  parameter  of  particular  im- 
portance to  an  absolute  determination  of  flux   is  the  diffusion  coefficient 
D;     for  as  can  be  seen  in  Eqs .    20,    24,    and  25,    it  serves  to  scale  the 
value  of  the  absolute  flux.     A  dilemma  arises  when  considering  the  best 
value  to  use  for  the  diffusion  coefficient  of  graphite  as  calculated  values 
and  experimental  values  may  differ  as  much  as  10  per  cent  (see  Table  7). 

D    can  be  calculated,   in  principle,    from  expressions  such  as  Eqs. 
37,    38,    40  and  41  below  (15). 


D  ,    1 j, (37) 


In  this  expression,     |j.      is  the  average  cosine  of  the  scattering  angle 
in  the  laboratory  system,    given  by  2 /3A    where    A    is  the  atomic 
weight. 

For  weak  absorbers  where    S    /S    «  1,   this  equation  reduces  to 

a'      s  m 


(38) 


3S.(»  -  ?„    ) 
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Defining  the  transport  mean  free  path  as 

\r     =  '       -  '  <39) 

the  diffusion  coefficient  becomes 

D    =    -*-  (40) 

From  a  development  of  the  significance  of  the  diffusion  length 
for  monoenergetic  neutrons  in  a  weak  absorber  (15) 

D    =    2a  L2  ,  (41) 

where    Z       is  evaluated  at  the  mean  energy  of  the  thermal 

spectrum  (43). 

An  often  quoted  value  of    D    is  calculated  from  Eq.    40  where 

X        is  calculated  from  Eq.    39.       X       appears  to  be  known  with 
tr  s 

some  certainty  as  0.385  cm"      (21);     and    (1  -  }I    )    is  calculated 

as  0.  944,(12).      This  results  in  a    X        in  the  range  from  2.71  cm 

to  2.75  cm  and  corresponding  diffusion  coefficients  in  the  range 

from  0.  903  cm  to  0.917  cm,    respectively.     This  has  been  the  range 

of  D    values  used  in  previous  pile  standardizations  (II,    35). 

It  has  not  been  feasible  to  calculate  adequately    D    from  Eq.    41; 

for  until  recently,    accurate  direct  measurements  of    S       were  dif- 

a 

ficult  to  make. 

Direct  measurements  of    D    thus     X        made  since   1955,    (6,    29), 
by  pulsed  source  techniques  disagree  somewhat  with  the  calculated 
values  outlined  above.      These  two  independent  measurements  of  X 
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have  been  reported  as   2.  55  ±  0.  09  and  2.  58  i  0.  02,   both  disagreeing 

considerably  with  the  values  for    X^       calculated  above. 

'  tr 

Amaldi  (2)  claims  that  one  cannot  rely  on  the  values  of  the  dif- 
fusion parameters  computed  from  scattering  and  absorption  cross 
sections  of  the  atoms  of  the  medium.     He  points  out  that  th'ese  various 
relationships  for  the  diffusion  parameters  are  useful  however,    since 
they  provide  the  possibilities  of  checking  in  detail  the  fulfillment  of 
the  assumptions  underlying  their  derivation. 

From  transmission  experiments,    one  measures     £       not    X     . 

X^       is  found  from    X      by  involvement  of    u    ,    the  calculation  of 

tr  s       '  'o 

which  depends  upon  the  angular  distribution  of  the  scattered  neutrons 

and  this  depends  in  a  rather  complicated  way  upon  the  molecular 

structure  of  the  medium. 

Singwi  and  Kothari  (36),    in  pointing  out  this  discrepancy  between 

the  calculated    X,       a:id  the  measured    Xt      ,    offer  no  conclusion  as  to 
tr  tr 

the  best  value  to  use.      They  say  that  the  theory  in  the  case  of  graphite 
is  still  in  an  unsatisfactory  state  and  that  a  better  theoretical  calcu- 
lation must  await  a  better  model  for  the  vibrational  spectrum  of 
graphite.     They  suggest  that  a  direct  measurement  of  the  angular  dis- 
tribution of  the  scattered  neutrons  in  the  thermal  energy  group  would 
be  valuable . 

On  the  basis  of  this  review  and  from  communication  with  the  Critical 
Mass  Physics  Section  of  Hanford  Laboratories  (9),    it  was  felt  that  the 
most  reasonable  value  of  the  diffusion  coefficient  would  be  found  from 

the  directly  measured  values  of    X       .     Thus     D    was  taken  as 

tr 

2.58/3  a  0.860  cm    at  a  graphite  density  of  1.60  gm/cm3. 
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A  paper  by  P.    F.   Nichols  (26)  gives  a  direct  measurement  of  the 
graphite  absorption  cross  section  in  the  Physical  Constants  Testing 
Reactor.     The  reactor  grade  graphite  (GAF)  tested  had  a  2200  m/sec 
cross  section  of    tr      =  3.80  ±  0.  04  mb    including  all  impurities. 
Using  this  and  the  diffusion  length  measured  in  the  KSU  pile  of  55.  3  ± 
0.3  cm,    a  diffusion  coefficient  of  0.828  cm  is  obtained  from  Eq.    41, 
thus  further  indicating  that  the  diffusion  coefficient  obtained  from 
Eq.    38  may  be  an  overestimate. 

Table  7  lists  a  summary  of  the  diffusion  parameters  of  graphite. 

Age.       In  the  finite  pile  flux  calculations,   the  Fermi  age  measured 

by  Steichen  (39)  has  been  used,    whereas  in  the  standard  pile  absolute 

flux  calculations,    a  combination  of  Gaussian  ranges  has  been  used 

instead  of  a  Fermi  age.     In  either  case,    to  evaluate  the  source  for  the 

thermal  diffusion  equation,   the  Fermi  age  or  a  Gaussian  range  to 

indium  resonance  must  be  extended  to  represent  the  age  or  range  at 

thermal  energy  by  Eq.    17  which  resulted  from  Fermi's  slowing  down 

theory 

Eln 

^th  s 

where 

rf,     .     =    4t  (E    .,    E ..    ).  (43) 

th-i  '    oi         th  '  x      ' 

Using  the  following  diffusion  parameters  appropriate  to  this 
work, 
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D     =     0. 

.860  cm 

E 

In 

=     1 . 44  ev 

2     =    0, 

s 

385  cm"' 

Eth 

=    0.025  ev 

i  =  o. 

158 

the  age  integral  above  has  a  value  of  57.  3  cm  . 

The  Reactor  Handbook  (30)  gives  an  estimated  age  from    indium 
resonance  energy  to  thermal  energy  of  53  cm2.       Weinberg  and 
Wigner  (43)  estimate  the  age  from  resonance  to  thermal  as  39  cm 
by  means  of 


~In 
(KJ  =  r  (ET     )  +  — (44) 


Lth       3  E 


J  J 

A  study  of  the  slowing  down  of  neutrons  in  graphite  by  Duggal 

and  Martelly  (13)  of  France  states  that  the  most  probable  value  of 

the  effective  slowing  down  age  between  indium  resonance  and  thermal 

equilibrium  is  49  cm  .     This  value  of  49  cm    was  used  for  this  work 

since  it  seemed  to  lie  in  between  the  extreme  values  and  it  agreed 

quite  well  with  that  given  by  the  Reactor  Handbook  (30). 

Thermalization  theory  (43)  was  not  applicable  in  this  case  since 
there  were  no  heavy  absorbers  such  as  uranium  to  harden  the  spectrum. 

The  age  of  neutrons  from  a  Pu-Be  neutron  source  to  indium  reson- 
ance energy  as  measured  by  Steichen  (39)  was  taken  as  419  cm    at  a 
graphite  density  of  1 .  6  gm/cm1.      This  number  was  the   result  of  an 
experimental  determination  of  422  and  a  theoretical  calculation  of  416. 
Extrapolating  this  age  to  thermal  energy  with  the  estimated  age  dif- 
ference of  49  cm    gave  the  age  of  neutrons  from  a  Pu-Be  source  to 
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thermal  energy  as  468  cm2.     At  the  effective  density  of  the  KSU  pile 
of  1.678  gm/cm3,   the  age  was  381  cm2  to  indium  resonance  and 
426  cm2  to  thermal. 

To  extend  any  indium  resonance  Gaussian  range  to  the  thermal 
Gaussian  range,    the  expression 

I«£    =    I]*±    ♦  (49)2  (45) 

was  used. 

Diffusion  Length.        The  diffusion  length  was  experimentally 
determined  to  be  52.76  ±0.30  cm    by  Kaiser  (22)  at  the  effective  density 
of  the  KSU  pile.     At  a  density  of  1 .  6  gm/cm3  this  would  be  55.  33  ±  0.  30 
cm.     This  number  is  the  result  of  the  average  of  ten  determinations  of 
diffusion  length,    each  determination  getting  1,  000  counts  with  the  BF3 
probe  in  each  of  13  selected  pile  positions.     The  model  from  which  the 
diffusion  length  was  determined  assumed  a  thermal  point  source. 

It  would  be  a  difficult  task  to  propagate  the  errors  of  each  individual 
parameter  to  give  the  confidence  limits  in  the  standard  pile  flux  as 
determined  by  Eq.    20.     Considering  the  uncertainty  in  the  diffusion 
coefficient  to  be  ±  5  per  cent,    (  estimated  from  the  difference  between 
0.860  and  0.910),   the  uncertainty  in  the  thermal  flux,    considering  only 
the  effect  of  the  source  strength  and  the  diffusion  coefficient,    amounts 
to  ±  12  per  cent.     Error  in  the  diffusion  length  (55.  33  ±  0.  30),    or 
estimated  error  in  the  age  (39)  or  a  Gaussian  range  should  have  a  small 
effect  on  the  flux  compared  with  the  effect  of  the  error  in  the  source 
strength  or  the  diffusion  coefficient.     A  study  of  the  variation  of  standard 
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flux  with  variation  of  diffusion  length  is  presented  in  Fig.    li 


Intercalibration  with  Argonne  Standard  Pile 
and  Hanford  Standard  Pile 


Table  8  presents  standardization  constants  of  KSU  gold  foils  as 
determined  from  the  irradiation  of  the  bare  foils  in  known  flux  posi- 
tions of  the  Argonne  Standard  Pile  and  the  Hanford  Standard  Pile. 
The  foils  were  irradiated  for  approximately  two  half-lives  and  then 
sent  to  Kansas  State  University  where  the  saturation  activity  of  the 
foils  was  determined  by  use  of  the  KSU  standard  foil  counting  system 
described  in  the  section  on  experimental  facilities.     On  the  basis  of 
the  saturation  activity  of  the  foil  and  the  knowledge  of  the  flux  to 
which  the  bare  foil  was  exposed,    a  standardization  constant  was  deter- 
mined for  the  foil  according  to 

(nv)th 

■Sh<«0-x-r  (46) 

sat 
The  results  of  the  subsequent  use  of  these  foils  in  the  KSU  Standard 
Pile  may  be  seen  in  Table  9. 

It  must  be  pointed  out  that  the  pile  at  Hanford  was  standardized 
with  a  diffusion  coefficient  of  0.  910  cm  ( 1 1)  and  the  pile  at  Argonne 
was  standardized  with  a  diffusion  coefficient  of  0.  903,    (32,    35). 
The  KSU  pile  was  standardized  on  the  basis  of  a  diffusion  coefficient 
of  0.860  cm  since,    as  previously  pointed  out,   this  was  thought  to  be 
a  better  estimate  of  the  diffusion  coefficient  for  graphite.     To  inter- 
calibrate  the  fluxes  of  the  various  standard  piles  on  the  basis  of  a 
consistent  diffusion  coefficient,    the  Argonne  flux  was  multiplied  by 
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Fig.  18.    Effect  of  diffusion   length  variation  on    (<SU    Standard    Pile 
flux. 
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0.  903/0.860  and  the  Hanford  flux  was  multiplied  by  0.  910/0860. 
Such  an  adjustment  led  to  the  flux  listed  in  Column  6  of  Table  9 
which  was  the  flux  to  be  compared  with  the  KSU  assigned  value 
of  the  absolute  flux. 

One  cause  of  the  unsatisfying  intercomparison  was  the  assump- 
tion expressed  in  Eq.    46  that  the  thermal  flux  is  proportional  to  the 
saturation  activity  of  a  bare  gold  foil.     Since  gold  has  its  peak  cross 
section  at  an  epi thermal  energy,    a  considerable  portion  of  the  acti- 
vation of  the  foil  near  the  source  was  due  to  epi  thermal  neutrons. 
The  epithermal  content  of  the  fluxes  in  the  intercalibrated  cases  was 
not  necessarily  the  same  since  a  KSU  neutron  source  has  a  different 
energy  spectrum  and  since  the  irradiations  took  place  at  different 
distances  from  the  source.     The  correct  procedure  for  the  intercali- 
bration  of  standard  pile  fluxes  would  be  the  exchange  of  both  bare  and 
cadmium  covered  gold  foils  since  the  thermal  flux  is  more  correctly 
proportional  to  the  cadmium  difference  saturation  activity. 

The  most  reliable  intercomparisons  made  with  this  bare  foil 
exchange  are  those  made  at  the  positions  fartherest  from  the  source, 
even  though  at  these  points  there  is  still  an  epithermal  flux  contri- 
bution. 

Theoretical  Standardization  of  BF3  Neutron  Probe 

A  theoretical  standardization  constant  for  the  BF3  neutron  detector 
was  determined  according  to  a  development  by  Price  (28).     Based  on 
the  manufacturer's   specifications  for  the  BF3  probe  described  in  the 


70 


section  on  experimental  facilities,    the  theoretical  standardization  con- 
stant of  the  probe  was  calculated  to  be  0.872  neutron  flux  per  count 
per  minute.     The  standardization  constant  determined  from  this  thesis 
was   1.  1578.     Considering  the  average  of  1 .  1578  and  0.  872  to  be  cor- 
rect,   this  amounts  to  a  discrepancy  between  the  two  of  28  per  cent  with 
the  theoretical  standardization  constant  giving  the  lower  flux. 

All  assumptions  made  in  the  theoretical  calculation  appear  to  be 
reasonably  correct.     However,    experience  at  Kansas  State  Univer- 
sity with  the  construction  of  BF3  neutron  detectors  (25)  indicates  that 
such  important  specifications  as  the  gas  pressure  and  the  active  volume 
are  difficult  to  determine  accurately.     Perhaps  the  specifications  given 
by  the  manufacturer  were  meant  to  be  nominal. 

Scintillation  Probe  Standardization 

The  standard  pile  analysis  gave  a  standardization  constant  of  0.8465 
neutron  flux  per  count  per  minute  for  the  scintillation  probe.      The  in- 
strument manufacturer  assigned  a  sensitivity  of  0.  014  counts  per  second 
per  unit  neutron  flux  to  the  probe.     Expressing  the  standardization  con- 
stant of  0.8465  in  terms  of  counts  per  second  per  unit  neutron  flux 
results  in  a  sensitivity  of  0.0197.     Considering  the  average  of  0.0197 
and  0.  014  to  be  correct,    a  34  per  cent  discrepancy  appears  to  be  involved 
with  the  manufacturer's  assigned  sensitivity  giving  the  higher  flux. 

It  must  be  pointed  out  that  the  manufacturer's  assignment  of  the 
sensitivity  was  based  on  a  calibration  with  the  Argonne  Standard  Pile. 
In  order  to  compare  the  KSU  standard  flux  with  the  flux  indication  from 
an  Argonne  standardized  detector,    the  detector  indication  must  be  ad- 
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justed  in  order  to  base  the  comparison  on  consistent  diffusion  coeffi- 
cients.    Such  a  consideration  would  decrease  the  manufacturer's  as- 
signed value  of  the  sensitivity  from  0.  014  to  0.  0133,    giving  a  larger 
39  per  cent  discrepancy  with  the  KSU  assigned  sensitivity  value  of 
0.  0197. 

Inspection  of  Table  6  or  Fig.    17  reveals  that  tl       scintillation  probe 
gave  the  least  reliable  standardization  constant  of  the  various  detectors. 
The  lower  standardization  constants  near  the  source  suggest  that  a 
standardization  using  a  cadmium  difference  technique  would  give  more 
reliable  results  for  the  scintillation  probe. 


The  Use  of  Gaussian  Ranges  and  the  Fermi 
Age  to  Describe  Slowing  Down 


Steichen  (39)  has  shown  that  the  age  from  a  polyenergetic  neutron 


source  may  be  theoretically  represented  by 
N 


*Sv, 


(47) 


i  =  l 

where    f.     represents  the  fractional  contribution  of    r. ,     the  age 

from  the    i         energy  group  of  the  source  to  an  energy  such  as 

indium  resonance  energy.     Writing    T.     in  terms  of    r.     according 

to  the  defining  equation  of    r. ,     Eq.    47  may  be  written  as 
N 


S* -3- 


(48) 


i  =  l 
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This  value  of  the  age,   hereafter  referred  to  as  the  equivalent  age, 
has  been  calculated  for  each  of  the  empirical  fits  and  tabulated  in 
Column  16  of  Table  5.     Inspection  of  this  column  reveals  that  even 
though  the  equivalent  ages  are  near  381  cm2.  Steichen's  (39)  value 
of  the  age  to  indium  resonance,    they  could  not  be  used  to  determine 
accurately  the  age.    Steichen  shows  that  the  value  of  the  Fermi  age 
is  determined  to  a  great  degree  by  the  extrapolation  to  infinity  of  the 
measured  resonance  distribution  beyond  the  last  measured  point  at 
97.  12  cm.      This  variation  of  the  equivalent  ages  with  choice  of  em- 
pirical fit  may  be  a  result  of  the  different  extrapolations  beyond 
97.  1  2  cm  given  by  various  empirical  fits  . 

It  was  of  interest  to  study  the  variation  of  the  calculated  standard 
pile  flux  that  resulted  from  different  choices  of  the  indium  resonance 
data  fit  given  in  Table  5.     Two  different  fits  were  used  to  calculate  the 
standard  pile  flux.     One  was  the  appointed  best  fit  which  had  an  equi- 
valent age  of  411.  9  cm2.     The  other  was  a  three  term  fit  of  the  indium 
resonance  data  corrected  for  geometry  and  foil  weights,   but  not  cor- 
rected for  neutron  absorption  by  higher  energy  indium  resonances. 
This  latter  fit  had  an  equivalent  age  of  378.  1  cm2.     Even  though  the 
equivalent  ages  and  the  individual  terms  of  the  two  fits  were  different, 
the  small  difference  in  the  calculated  fluxes  was  difficult  to  show  clear- 
ly in  a  plot  of  the  two  fluxes.      The  maximum  difference  was  a  little 
greater  than  2  per  cent. 

Use  of  the  three  Gaussian  ranges  allows  a  more  accurate  descrip- 
tion of  the  slowing  down  density  than  does  the  use  of  the  Fermi  age. 
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As  can  be  seen  from  the  form  of  Eq.    8,    use  oi  the  Fermi  age  to  describe 
the  slowing  down  distribution  would  yield  a  straight  line  on  Fig.    14.     The 
use  of  three  Gaussian  ranges,    hovvi  ,  er,    provides  an  adequate  descrip- 
tion of  the  actual  measured  slowing  down  distribution. 

Further,    it  must  be  stated  that  the  use  of  the  three  Gaussian  ranges 
to  determine  the  slowing  down  to  thermal  energy  will  allow  a  more 
accurate  description  of  the  thermal  neutron  flux  than  will  the  use  of 
the  Fermi  Age.      This  conclusion  has  not  been  explicitly  demonstrated 
but  it  follows  from  the  fact  that  a  more  accurate  description  of  the 
slowing  down  of  neutrons  to  thermal  energy  will  allow  a  more  accurate 
source  boundary  condition  to  the  therinal  diffusion  equation. 

Study  of  Source  Boundary  Condition  in  Finite  Pile 

Table   10  is  a  summary  of  flux  calculations  in  the  study  of  the  source 
boundary  condition  to  the  thermal  diffusion  equation  in  the  finite  pile 
geometry.     Column  2  lists  the  flux  using  fast  point  source  theory  as 
obtained  from  the  computer  code  described  in  Appendix  E  which  numer- 
ically evaluates  Eq.    24.     Column  3  lists  the  flux  using  thermal  point 
source  theory  as  obtained  from  the  computer  code  described  in  Appen- 
dix F    which  numerically  evaluates  Eq.    25.     Column  4  lists  the  exper- 
imental absolute  flux  determined  by  the  standardized  BF3  probe.     All 
measurements  and  calculations  were  made  with  a  source  having    a 
strength  oi    I    71x10    neutrons  per  second. 

Figures  20  are  isometric  plots  of  the  flux  values  in  Table   10. 
The  absolute  flux  is  plotted  along  the  vertical  axis  as  a  function  of 


Table  1C.  Study  of  the  flux  from  two  assumed  source  boundary  conditions  for 
the  thermal  diffusion  equation  in  a  finite  pile.   Q0  =  1.71  x  1  06  n/sec. 
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Distance  from 

Thermal  flux 

Thermal  flux 

Thermal  flux 

source 

from  fast 

from  thermal 

from  standardized 

X 

:     Y     :     Z 

point  source 

point  source 

BF3  probe 

inches 

n/cmz-sec 

n/cm  -sec 

n/cm2-  sec 

0 

0          0 

2,238.0 

4 

2,162.0 

13,150.0 

2,432.0 

8 

1,951.0 

5,356.0 

2,190.0 

12 

1,647.0 

2,879.0 

1,794.0 

16 

1,304.0 

1,716.0 

1,423.0 

20 

966.0 

1,061.0 

1,115.0 

24 

659.4 

650.5 

666.9 

28 

391.6 

363.2 

389.8 

32 

154.3 

139.1 

167.6 

0 

8          0 

1,951.0 

5,356.0 

2,094.0 

4 

1,886.0 

4,594.0 

2,037.0 

8 

1,706.0 

3.177.0 

1,820.0 

12 

1,445.0 

2,101.0 

1,508.0 

16 

1,150.0 

1,377.0 

1,115.0 

20 

855.4 

896.5 

842.6 

24 

586.8 

569.6 

585.6 

28 

349.9 

320." 

338.7 

32 

138.2 

123.8 

136.4 

0 

16          0 

1,304.0 

1,716.0 

1,335.0 

4 

1,263.0 

1,615.0 

1,293.0 

8 

1,150.0 

1,373.0 

1,159.0 

12 

982.7 

1,078.0 

982.4 

16 

790.8 

799.9 

771.7 

20 

596.7 

566.9 

578.0 

24 

4H.5 

378.3 

404.4 

28 

249.8 

222.6 

257.0 

32 

99.29 

87.41 

100.3 

0 

24          0 

659.4 

650.5 

621.7 

4 

640.5 

629.5 

511.3 

8 

586.8 

565.8 

566.5 

12 

507.7 

474.6 

482.2 

16 

414.5 

378.5 

396.0 

20 

317.7 

283.0 

300.7 

24 

224.2 

196.5 

219.0 

28 

136.8 

118.2 

132.5 

32 

54.74 

47.33 

54.74 

0 

32          0 

154.3 

126.3 

147.1 

4 

150.1 

134.8 

142.8 

8 

138.2 

125.9 

133.4 

12 

120.5 

106.2 

113.4 

16 

99.29 

86.61 

94.42 

20 

76.87 

67.36 

68.31 

24 

54175 

47.25 

52.33 

28 

33.64 

28.60 

32.53 

32 

13.53 

11.58 

14.39 
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Table 

10   (cont.) 

Distance  from 

Thermal  flux 

Thermal  flux 

Thermal  flux 

source 

from  fast 

from  thermal 

from  standardized 

X 

:     I     :     Z 

point  source 

point  source 

BF3   probe 

inches 

n/cm2-sec 

n/cm  -sec 

n/cm  -sec 

0 

0          8 

1,951.0 

5,356.0 

2,148.0 

4 

1,886.0 

4,594.0 

2,062.0 

8 

1,706.0 

3,177.0 

1,815.0 

12 

1,445.0 

2,101.0 

1,516.0 

16 

1,150.0 

1,377.0 

1,188.0 

20 

855.4 

896.5 

861.4 

24 

586.8 

564.6 

605.5 

23 

349.9 

320.8 

339.2 

32 

138.2 

123.8 

H8.2 

0 

8          8 

1,706.0 

3,177.0 

1,775.0 

4 

1,650.0 

2,889.0 

1,721.0 

8 

1,495.0 

2,258.0 

1,536.0 

12 

1,271.0 

1,623.0 

1,266.0 

16 

1,015.0 

1,130.0 

1,012.0 

20 

759.3 

764.0 

747.4 

24 

523.3 

492.5 

514.9 

28 

313.2 

283.8 

318.7 

32 

124.0 

110.3 

126.4 

0 

16          8 

1,149.0 

1,373.0 

1,162.0 

4 

1,114*0 

1,312.0 

1,105.0 

s 

1,015.0 

1,127.0 

1,002.0 

12 

870.7 

909.9 

896.7 

16 

703.3 

689.3 

692.9 

20 

532.8 

496.3 

528.0 

24 

371.8 

334.7 

371.9 

28 

224.8 

198.5 

224.0 

32 

89.52 

78.25 

91.62 

0 

24          8 

586.8 

565.8 

564.4 

4 

570.2 

541.0 

542.4 

8 

523.3 

491.3 

498.0 

12 

453.9 

417.4 

434.2 

16 

371.8 

335.2 

362.4 

20 

285.9 

253.3 

273.5 

24 

202.4 

176.6 

192.6 

28 

123.8 

106.9 

122.5 

32 

49.63 

42.69 

49.94 

0 

32          8 

138.2 

125.9 

128.86 

4 

134.5 

121.8 

128.32 

8 

124.0 

111.3 

119.75 

12 

108.3 

96.36 

103.08 

16 

89.52 

78.03 

84.44 

20 

69.50 

60.43 

66.96 

24 

49.63 

42.84 

47.53 

28 

30.56 

26.19 

29.99 

32 

12.30 

10.52 

13.89 
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Table  10  (cont.) 

Distance  from 

Thermal  flux 

Thermal  flux 

Thermal  flux 

source 

from  fast 

from  thermal 

from  standardized 

X     :     I     :     Z 

point  source 

point  source 

BF3  probe 

inches 

n/cm2-sec 

n/cm2-sec 

n/cm  -sec 

0          0         16 

1,304.0 

1,716.0 

1,330.0 

4 

1,263.0 

1,615.0 

1,269.0 

8 

1,150.0 

1,373.0 

1,167.0 

12 

982.7 

1,078.0 

992.5 

16 

790.8 

799.9 

774.3 

20 

596.5 

566.9 

594.4 

24 

419.5 

378.3 

412.4 

28 

249.8 

222.6 

252.4 

32 

99.29 

87.41 

103.7 

0          8        16 

1,149.0 

1,373.0 

1,118.0 

4 

1,114.0 

1,312.0 

1,092.0 

8 

1,015.0 

1,127.0 

1,010.0 

12 

870.7 

909.9 

870.7 

16 

703.3 

689.3 

696.7 

20 

53.28 

496.3 

529.3 

24 

37.18 

334.7 

365.6 

28 

224.8 

198.5 

224.30 

32 

89.52 

78.25 

92.43 

0        16        16 

790.8 

799.9 

753.4 

4 

767.9 

772.3 

767.9 

8 

703.3 

690.2 

686.9 

12 

608.2 

576.8 

597.4 

16 

496.2 

457.8 

482.5 

20 

480.1 

341.2 

367.9 

24 

268.0 

235.7 

274.1 

28 

163.4 

141.7 

166.8 

32 

65.41 

56.37 

66.69 

0        24        16 

414.5 

378.5 

398.3 

4 

403.4 

365.8 

335.9 

8 

371.8 

334.4 

■359.4 

12 

324.6 

290.7 

310.7 

16 

268.0 

234.8 

255.3 

20 

208.0 

181.3 

198.9 

24 

U8.4 

128.6 

143.20 

23 

91.34 

78.90 

91.85 

32 

36.77 

31.48 

39.02 

0         32        16 

99.26 

86.61 

94.94 

4 

96.75 

84.95 

90.52 

8 

89.52 

74.81 

87.22 

12 

78.64 

66.14 

77.32 

16 

65.82 

57.40 

63.48 

20 

51.15 

43.31 

49.32 

84 

36.77 

31.66 

36.41 

-  > 

22.76 

19.82 

22.15 

! 

9.02 

7.896 

10.74 

Table  10  (coat.) 


Distance  from 

source 
X     :     I     :     Z 

Thermal  flux 

from  fast 
point  source 

Thermal  flux 
from  thermal 
point  source 

Thermal  flux 
from  standardized 
BF3    probs 

inches 

n/cm2-sec 

n/cm2-sec 

n/cm2-sec 

0 

4 

8 
12 
16 
20 
24 
28 
32 

0 
4 
8 

12 
16 
20 
24 
28 
32 

0 
4 
8 
12 
16 
20 
24 
28 
32 

0 
4 
8 
12 
16 
20 
24 
23 
32 

0 

4 

8 

12 

16 
20 

;  % 

28 
32 


24 


8   24 


16 


24 


24   24 


32   24 


659.4 

650.5 

650.0 

640.5 

679.5 

630.2 

586.8 

565.8 

578.4 

507.7 

474.6 

502.3 

4H.  5 

378.5 

413.6 

317.7 

283.0 

317.7 

224.2 

196.5 

227.6 

136.8 

118.2 

H2. 6 

54.74 

47.33 

61.80 

586.8 

565.8 

554.1 

570.2 

541.0 

553.7 

523.3 

491.3 

517.3 

453.9 

417.4 

445.5 

371.8 

335.2 

363.6 

285.9 

253.3 

277.9 

202.4 

176.6 

197.1 

129.3 

106.9 

120.6 

49.63 

42.69 

53.29 

414.5 

378.5 

405.9 

403.4 

365.8 

386.4 

371.8 

334.4 

349.7 

324.5 

290.7 

313.5 

268.0 

234.8 

265.4 

208.0 

181.3 

205.5 

148.4 

128.6 

H4.5 

91.34 

78.90 

93.53 

36.77 

31.48 

38.71 

224.2 

196.5 

215.2 

218.5 

189.3 

207.7 

202.4 

175.3 

196.7 

178.1 

154.6 

168.3 

US.  3 

128.0 

144.3 

116.3 

100.2 

110.2 

83.75 

71.77 

82.53 

51.92 

44.84 

54.10 

20.99 

18.18 

24.02 

54.75 

47.25 

52.60 

53.42 

45.21 

51.27 

49.63 

41.95 

47.04 

43.86 

37.11 

44.00 

36.77 

32.30 

37.11 

28.93 

24.25 

29.23 

20.99 

17.65 

21.40 

13.07 

11.53 

14.33 

5.30 

4.67 

7.56 

Table  10  (concl.) 
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Distance  from 

Thermal  flux 

Thermal  flux 

Thermal  flux 

source 

from  fast 

from  thermal 

from  standardized 

X 

:     2     : 

Z 

point  source 

point  source 

BF3   probe 

inches 

n/rmJ-sec 

n/cm  -sec 

n/cm2- sec 

0 

0 

32 

154-3 

139.1 

149.4 

4 

150.1 

134.8 

149.1 

8 

138.2 

125.9 

138.1 

12 

120.5 

106.2 

121.9 

16 

99.29 

86.61 

95.63 

20 

76.87 

67.36 

74-66 

24 

54.75 

47.25 

54. 42 

23 

33.64 

28.60 

34.13 

32 

13.53 

11.53 

14.47 

0 

3 

32 

134.9 

125.9 

126.06 

4 

134.5 

121.8 

128.28 

8 

124.0 

111.3 

122.73 

12 

108.3 

96.36 

100.8 

16 

89.52 

78.03 

83.99 

20 

69.50 

60.43 

68.66 

24 

49.63 

42.84 

51.52 

23 

30.56 

26.19 

28.13 

32 

12.30 

10.52 

13.69 

0 

16 

32 

99.29 

86.61 

99.10 

4 

96.75 

84.95 

88.90 

8 

89.52 

74.81 

84.50 

12 

78.64 

66.14 

76.80 

16 

65.42 

57.40 

65.57 

20 

51.15 

43.31 

50.51 

24 

36.77 

31.66 

35.95 

28 

22.76 

19.82 

23.62 

32 

■ 

9.195 

7.90 

10.08 

0 

24 

32 

54.75 

47.25 

52.74 

4 

53.42 

45.21 

47.91 

8 

49.63 

41.95 

47.27 

12 

43.86 

37.11 

40.81 

16 

36.77 

32.30 

35.31 

20 

28.98 

24.25 

27.73 

24 

20.99 

17.65 

20.64 

28 

13.07 

11.53 

12.09 

32 

5.297 

4.67 

5.89 

0 

32 

32 

13.53 

11.58 

12.39 

4 

13.21 

12.64 

8 

12.31 

10.53 

13.35 

12 

10.92 

11.61 

16 

9.195 

7.90 

8.949 

20 

7.279 

8.660 

24 

5.297 

4.67 

5.708 

28 

3.304 

4.006 

32 

1.343 

1.22 

1.424 
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Fig.  19.    Thermal   flux  study  in   finite  pilo. 
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x,   y  and  z  with  each  figure  representing  a  different    z  level  in  the 
finite  pile.     Figure   19  is  a  semi -logarithmic  plot  of  the  flux  as  a 
function  of    x    alone,      z  and    y    being  zero.     Figure   19  clearly  shows 
the  difference  between  the  flux  from  fast  source  theory  and  the  flux 
from  thermal  source  theory.     The  fast  source  calculations  used  the 
single  value  of  the  age  of  neutrons  from  a  Pu-Be  source  to  thermal 
energy.     Figure  19  shows  the  effect  of  using  the  three  range  groups, 
or  an  analytical  source  expression,   determined  from  the  standard 
pile  analysis  instead  of  the  single  value  of  the  age.      This  curve,    cal- 
culated by  Eq.    24,    resulted  from  the   superposition  of  the  flux  con- 
tributions from  three  different  monoenergetic  sources,    each  source 
having  its  respective  Gaussian  range  and  fractional  contribution 
determined  empirically  from  Eq.    36  and  extended  to  thermal  energy 
by  Eq.    45.     It  can  be  seen  that  the  breaking  up  of  the  source  into  three 
energy  groups  by  the  use  of  Gaussian  ranges  does  give  an  improved 
effect. 

Studies  of  the  variation  of  the  flux  in  fast  source  theory  and  thermal 
source  theory  with  the  variation  of  various  parameters  are  given  in 
Figs.    21  through  23.     Graphs  showing  the  variation  of  flux  with  vari- 
ation of  source  strength  or  diffusion  coefficient  are  not  given,    since 
both  of  these  parameters  are  merely  scaling  coefficients  to  the  flux. 

The  values  of  the  source  strength  and  the  diffusion  coefficient  do 
not  affect  the  various  plots  in  Figs.    19  and  20  with  respect  to  each 
other  since  they  were  all  based  on  consistent  values. 
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Fig.  21.   Effect   of  variation  of  Fermi    age  on  thermal  flux    from 
fast    point   source    in  finite  pile. 
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Fig.  22.   Effect   of  variation  of    diffusion  length  on  thermal   flux 
from    fast  point  source    in     finite    pile. 
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Fig. 23.    Effect  of  variation   of  diffusion  length  on   thermal    flux 
from   thermal    point  source   in  finite  pile. 
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Suggestions  for  Further  Work 

The  IBM-650  program  which  calculates  the  thermal  flux  from  fast 
point  source  theory  in  finite  geometry  offers  interesting  possibilities 
for  research  with  superposition  of  fluxes.  Sources  may  be  placed  at 
various  positions  in  the  pile  and  the  absolute  flux  may  be  determined 
experimentally  with  the  standardized  detectors.  This  experimental 
absolute  flux  may  be  compared  with  the  superimposed  fluxes  deter- 
mined from  machine  calculations  of  the  flux  from  individual  sources. 

It  is  evident  that  there  is  still  uncertainty  in  defining  the  absolute 
flux.     An  independent  calculation  of  absolute  flux  could  be  made  with 
the  absolute  counting  of  beta  particles  from  indium  foils.     A  recent 
outline  for  such  work  is  References  (17)  and  (23). 

An  approximate  method  for  an  absolute  flux  determination  using 
a  water  medium  has  been  outlined  by  Reference   1.      The  method  assumes: 

(a)  The  count  rate  of  a  thermal  neutron  detector 

such  as  the  BFj  probe,   is  directly  proportional 
to  the  thermal  flux  at  all  points  in  the  water 
medium. 

(b)  The  water  medium  is  infinite  and  spherical 
geometry  can  be  used.     Reference  (30)  con- 
siders the  water  volume  necessary  to  approxi- 
mate infinity. 

(c)  Neutron  capture  occurs  only  at  thermal  energies. 

(d)  2       for  water  is  accurately  known. 

(e)  Q0    is  accurately  known. 

Since  leakage  from  the  "infinite  medium"  is  negligible,    the  only  loss  of 
neutrons  is  due  to  capture  at  thermal  energies.     It  can  then  be  written 
in  spherical  coordinates  that 
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o=    \  Sa4>(r)  4,rr2dr  (49) 


L'all  space 

From  the  assumption  that    $  =  KR(r)    where     r    is  the  count  rate, 
K  may  be  written  as 

Qo 

K     =    -a, (c'f; 


r«. 


2    R(r)  r2dr 


l<     may  be  experimentally  determined  and  the  integral  may  be 
evaluated  numerically. 

Some  work  may  be  done  in  attempting  to  tie  up  the  foil  standard- 
izations in  graphite  with  the  absolute  flux  in  other  media  such  as 
water  or  paraffin.     A  first  approximation,    ignoring  any  difference 
in  the  energy  of  the  thermal  neutron  spectrum  in  graphite  and  water 

for  example,   would  require  the  determination  of    K°    ,     the 

th 
standardization  constant  independent  of  perturbing  effects  introduced 

in  the  section  on  data  analysis.      This  is  necessary  since  many  of 

the  perturbations  depend  upon  the  moderating  properties  of  the  medium. 

It  is  felt  by  this  author  that  a  complete  systematic  study  of  the 
foil  correction  factors  would  be  desirable.     In  addition  to  the  contribu- 
tion of  helping  to  make  possible  the  use  of  the  standardized  foils  in 
media  other  than  graphil    ,    this  study  would  be  necessary  in  determin- 
ing the  absolute  value  of  an  unperturbed  flux  from  absolute  beta  count- 
ing of  foils . 

An  interesting  addition  to  this  work  might  be  the  comparison  of 
the  standard  pile  flux  with  the  flux  obtained  from  the  finite  pile  solu- 
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tions  which  would  use  the  standard  pile  in  its  finite  dimensions.     It 
may  be  recalled  that  in  the  standard  pile  analysis,    the  pile  was 
assumed  to  have  an  infinite  height.     A  difficulty  of  interest  that  would 
arise  would  be  the  definition  of  the  boundary  condition  at  the  lower 
face  of  the  pile,   which  rests  on  concrete. 

In  the  discussion  of  the  results  it  was  stated  that  the  use  of 
Gaussian  ranges  to  determine  the  slowing  down  to  thermal  would  allow 
a  more  accurate  description  of  the  thermal  neutron  flux  than  would  the 
use  of  the  Fermi  age.     A  standard  pile  flux  calculation  using  Steichen's 
value  of  the  age  instead  of  the  three  Gaussian  ranges  would  show  this 
difference  in  thermal  flux. 
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APPENDICES 


APPENDIX   A 


Solutions  to  the  Fermi  Age  Equation  and 

the  Thermal  Diffusion  Equation  for  the 

Standard  Pile  Analysis 


The  geometry  for  the  standard  pile  solution  is  shown  in  Fig.    24. 


Fig.    24.     Geometry  for  the  standard  pile  solution. 

The  equations  solved  were  the  Fermi  age  equation 
9q 


V'q 


3t 


and  the  thermal  diffusion  equation 

DV24>     -     S   4>     +    S     =     0 

The  Fermi  age  equation  was  subjected  to  the  boundary  con- 
ditions that    q  =  0  at  x  =  ±a/2,       y  =  ±b/2,       and  q  —  0    as    |  zf  — 


(A-l! 


(A-2^ 


The  source  condition  was  that  all  source  neutrons  entered  the  system 
at    x  =  y  =  z  =  0,     at  an  age     t  =  0. 

This  condition  was  written  into  the  Fermi  age  equation  in  the  form 

V2q=|f  -Q„6(t)  6(x)6(y)  6(z).  (A-3) 

The  boundary  conditions  to  the  thermal  diffusion  equation  were 
that      (j>  =  0  at  x  =  ±a/2,       y  =  ±b/2,     and  cf>  —  0  as    |  z|   ->oo.     The 
source  condition  was  that    S  =  q(x,  y,  z,    T     ). 

Q0,     q  and   4>    were  expressed  as  double  Fourier  series  satis- 
fying the  boundary  conditions 

00  00 


m-1     n=  1 
and 


♦  -II    L>)-  ^  ">S^F  <a"4> 

m  =  l     n=l 

00  00 

»    =   I    Z        Vn  (Z'     COS  -1ILF-  COS  "Hp  (A"5) 

oo       oo 

Q.  6(t)  6(x)  6(y)  6(t)     =   £      J      Qmn  (*)  cos  JB^L  cos  Jill  (A-6) 

m  =  l     n  =  l 

where  m  and  n  are  summed  over  odd  integers.  Each  side  of 
Eq.  A-6  was  multiplied  by  orthogonal  functions  and  then  was  in- 
tegrated over  the  range  of  orthogonality  which  resulted  in 

4Q0 

Q0  5(t)  6(x)  6(y)  6(z)     =    6(t)  6(z)  cos  -5LE2L  cos  SgL  (A-7) 

ab 

Substitution  of  Eqs.   A-5  and  A-7  into  the  age  equation, 

Eq.   A-3  gave 
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CO  CO 

m-1  n~\ 


()q  o   Q  2    2  22 

+  — , +  )q  +  _^  6(t)  5(z) 


3r 


az< 


4Q0 
^mn    '    ab 


Taking  the  Fourier  transform  of  the  equation  for  a  particular 
m    and    n    and  letting 


(A-8) 


and 


led  to 


3      .    Si. 

n  b 


HJSL  +    («*         +   ^     +  «*  )   q     („)      =  "        6(T) 

„  m         r  n  /    i     \    /  aD         -v    / 

or 


(A-9) 
(A-10) 

(A-ll) 


where     q  («a)     represents  the  Fourier  transform  of    q  (z). 

The  solution  to  this  differential  equation  was 

4Q0      -(a1     +  (3*    +  M*  )r 
—  ,    >  m       r  n  ' 


(A-12) 


Taking  the  inverse  transform  of  Eq.    A-12  and  substituting 
it  into  Eq.    A-5  gave  the  slowing  down  density  along  the     z  axis  as 


-(a'     +  (1*    )t 
4Q0e        m  " 


**)*l      I 


7.*/4T 


m  =  l     n=l  ab  \l  4  ir  t 


(A-13) 


Substitution  of  Eqs.    A-4  and  A-5  into  the  thermal  diffusion 

E 
equation ,      Eq.    A  -  2,    and  designating    — =- —     as     k      gave 
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00  00 

1 1 

m  =  l     n  =  l 


8z« 


n  Tmn  D 


0. 


(A- 14) 


Taking  the  Fourier  transform  of  Eq.   A- 14  for  a  particular 

m    and    n    and  designating  the  Fourier  transforms  of  6  and  q 

■  Tmn  r 

by      <(>  (u)    and    q  (u),     respectively,    gave 


4-  M 


q  (") 


D(a2     +  /32    +  u2  ) 
*    m       'n  ' 


Taking  the  inverse  transform  of  Eq.    A-15  by  means  of  the 

convolution  theorem  (42),    after  Eq.    A- 12  had  been  used  to  express 

q  (w)  gave 

oo 


2Q0b         e 
mn 


m       r  m 


z"2/4t      -|z-zv|/b 


Dab  n/4ttt 


where 


V3 


°m+0_   +   KZ 


(A-15) 


dzv  (A-16) 


(A-17) 


Glasstone  and  Edlund  (15),    on  pages   185  to  187,    give  the 
development  showing  how  the  integral  in  Eq.    A  -16  may  be  expressed 
in  terms  of  the  error  function 

t 


erf(t) 


4V 


J     e"u2du 


(A-18) 


Writing  Eq.    A-16  in  terms  of  the  error  function  and  sub- 
stituting it  into  Eq.    A-4  gave  the  thermal  flux  along  the    z  axis  as 
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00  w  ^     T_ 

x,   ,         S  X  Q°bmn         *2 

*<z)    =   Z  Z  Dab  e 

jaa— _1  n  =  1 


s/b 


1+  erf 


n/t" 


2-v/T         2b 


z/b 


1   -  erf 


UsTT  2b 


(A- 19) 


If  difficulty  has  been  encountered  in  following  this  solution,    it 
may  be  helpful  to  study  References  (37)  and  (42). 
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APPENDIX   B 

Finite  Pile  Solution  to  the  Thermal  Diffusion 
Equation  Using  the  Distributed  Thermal 
Source  Found  from  the  Solution  to 
the  Fermi  Age  Equation. 


The  geometry  for  this  finite  pile  solution  is  shown  in  Fig.    25. 


Fig.    25.     Pile  geometry  for  fast  point  source  in  finite  pile. 


The  equations  solved  were    the  Fermi  age  equation 
3q 


V2q 


8  T 


and  the  thermal  diffusion  equation 

DV24>  -  S    +  +  S  =  0. 


(B-l) 


(B-2) 


The  Fermi  age  equation  was  subjected  to  the  boundary  con- 
ditions that    q  =  Oatx  =±a,       y  =  ±b,     and     z  =  ±  c.  The  source  condition 
was  that  all  source  neutrons  entered  the  system  at    x  =  x0,     y  =  y0 
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and    z  =  z0  ;     and  at  T  =  0. 

The  boundary  conditions  to  the  thermal  diffusion  equation  were 
that  <j>  =  0atx  =  ±a,  y  =  ±b,  and  z  =  ±c.  The  source  condition 
was  that    S  =  q(x,   y,  z,    T     ). 

By  assuming  that  the  variables  of    q    were  separable, 

q(x,   y,    z,    t)     =    M(x,   y,    z)N(t) 
the  age  equation  could  be  written 


(B-3) 


M  N      dT 


(B-4) 


Solutions  satisfying  the  boundary  conditions  and  the  two 

independent  differential  equations  were 
oo         oo        oo 


whe 


M(x,y,.)-2       I       1     A.^     Fm(x)  Gn(y)  Ho(z)  (B-5) 


m  =  l      n  =  l    o=l 


F (x)    =    cos  a      x 

•JST 


F     (x) 


sin  or      x 


(B-t) 


and 


00  00  00 


l(r)     =}      )       )       A" 
m  =  l    n=l    o=l 


B<  T 

mno 


(B-7) 


where 


B2  =  a2      + 

mno  m 


(B-8) 


104 


G   (v),    and  M   (z)    are  similar  to  F     (x)     except  that    b    and     8 
"  n   '  '  o     '  m     ' 

replace     a    and    a ,    respectively,    for  the    y   direction;     and    c     and 

r    are  used  for  the     z  direction. 

To  evaluate    A  in  the  resulting  expression  for    q. 

in  no  °        ■ 

00       co       oo  -B  T 

V"      s;-      V"  mno 

q     =     )  /A  F       (x)  G   (y)  H    (z)  e  (B-9) 

Lj      Lj      Lj       mno      m  *    '      nw'      o     ' 

m  =  l     n  =  l     o  =  l 
the  source  condition  that 

q(x,    y,    z,    0)  =     Q„    6[x  -  x„)  8(y   -  y0)  S(z   -   z0  )  (B-10) 

was  used.     Expressing  Eq.   B-10  with  the  use  of  Eq.    B-9,    multiplying 
both  sides  of  the  result  by  orthonorrnal  functions,    and  integrating 

over  the  range  of  orthonormality  led  to  the  expression  for    A 

'  mno 

Amno     "     Q»  Fm  (-0  )  Gn  (y0)     H^,  ).  (B-ll) 

The  expression  for  the  slowing  down  density  evaluated  at 

thermal  energy  was  then 

oo        co  oo  -B2        r 

V  \  \  mno    tli 

q(x,y,Z,   rA)=  2j  Z  L     Q(,Fm(x')Fm(xo)Gn(y)Gri(y0)H0(K)H(,(Zo)e 

m  =  I  n  =  l  o=l 

Assuming  that   cb    could  be  represented  by  the  same  ortho- 
normal  set  used  for    q,     a  solution  for  the  thermal  diffusion  was 


(B-1Z) 


<JQ  00  00 

f  y  y 


(x,   y,    z)  =     >         >         )  C  F     (x)  G    (y)  H     la).  (B-13) 

Lj      Lj      Lj  mno      m     '      nw '      o       '  l  ' 

m=l     n  =  l    o=l 
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Substitution  of  this  expression  and  the  expression  for  slowing 
down  density  evaluated  at  thermal  energy,    Eq.    B-12,    into  the  thermal 
diffusion  equation  yielded 


Qo  F     (x0)  G   (y„)  H     (z0)  "     mno     th 

C  = S2 - - e  (B-14) 

mno 

D(az      +  /32    +  72   +    k2  ) 
v    m         *n  o  * 

where 

k2   9    Sa  /D.  (B-15) 


The  thermal  flux  was  then 


V    X     Q0Fm(x)Fm(x0)Gn(y)Gn(y0)H„(z)H0(Zo)      "Bmno   "th 


'(X^2)  =     III 

1  D(  a2     +  II2    +    7*   +    K2) 

m  n  o 

(B-I6) 


m=l    n=l    o=l  D(  az     +  0Z    +   7Z   +    x' 

m  n  o 
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APPENDIX  C 


Finite  Pile  Solution  to  the  Thermal  Diffusion  Equation 
using  a  Point  Source  of  Thermal  Neutrons. 


The  solution  using  a  point  source  emitting  only  thermal  neutrons 
could  have  been  determined  by  merely  using    T       as  zero  in  Eq.    B-15. 
However,   it  was  found  that  such  an  expression  converged  quite  slowly. 
The  thermal  diffusion  equation  was  re-solved  using  the  geometry  shown 
in  Fig.    26. 


2c 


Fig.    26.     Pile  geometry  for  thermal  point  source  in  finite  pile. 

The  boundary  conditions  to  the  thermal  diffusion  equation  were 
that      <|>  =  0  at  x  ±a,     y  ±b,     and  z  =  0,      2c.       The  point  thermal 
neutron  source  was  considered  to  be  in  the    z  =  z0    plane  at 
x  =  x0,     and    y  =  y0.     The  source  was  introduced  through  the  boundary 
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conditions  that 


<t>i  (z<>)    =    <t>2  (zo) 


and 


-D 


d<(>i  d<j>2 

dz 


Qo 


where  $i  was  considered  to  be  the  flux  of  the  region  below 
the  z  =  z0  plane  and  <J>2  was  the  flux  of  the  region  above  the 
z0    plane. 

Assuming  that  the  variables  of    (J>    were  separable, 

4>  (x,   y,    z)     =    F(x)  G(y)  H(z) 

the  thermal  diffusion  equation  written  without  a  source  term. 


V2  <)>     -   K2    $     =    0 


could  be  written  as 

1      d2  F      ,      1      d2G 


dx! 


♦if?-*"* 

d  z 


A  solution  for  the  region  above  the     z0    plane  that  satis- 
fied the  boundary  conditions  and  the    diffusion  equation  was 


<t>(x,   y.  z)  =     )         ) 
m  =  l      n  =  l 


A         F     (x)  G   (y)  H       (z) 
mn      nV    '      nw '      mnv    ' 


where 


F     (x) 


1 


1 
"/a" 


cos  a      x 


for    m    odd 


for    m    even 


(C-l) 
(C-2) 


(C-3 


(C-4) 


(C-4) 


(C-6) 


(C-7) 
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and 


m  2a 


H         (z)    =     sin  h  7         (2c  -  z) 
mn  mn  ' 


(C-8) 


where 


72        =    a2-      +    |32   +    k2 
m         *n 


(C-9) 


G     (y)    is  similar  to    F     (x)    except  that    b,     n,    and    /3 
replace     a,   m,     and    a,    respectively. 

Expanding  the  source  term  in  a  similar  two-dimensional 
set  in    x    and    y, 


00  00 

Q0   6(x  -  x0  )  6(y  -  yo)     =    J,      I       Qmn  Fm(x)  G^ 


(y). 


m  =  l    n=l 


(CIO) 


multiplying  both  sides  by  the  orthonormal  functions  and  integrat- 
ing over  the  range  of  orthonormality  gave 


oo       oo 

50  6(x-  x„)  6(y  -  y„)    =    )        )        Q0  F       (x0)  G 

I— i       L—i  m  n 


(yo) 


(CXI) 


m  =  l     n  =  l 


Using  the  boundary  conditions,    Eq.    C-l  and  Eq.    C-2  with  Eq.   C-Il 
led  to  the  evaluation  of    A  for  the  region  above  the     z„    plane, 


Qq  Fm(*°)  Gn(yo)  sin  h   Tmn  (z0) 

D7  sinh     27  C 

mn  mn 


(C-12) 


The  expression  for  the  thermal  flux  in  the  region  above  the 
plane  of  the  source  was  then 

V       V       O  sinh  7        z0    sinh7       (2c-z 

^■Z>=Z      1     TPT    FmW  Fm(x0)  G(y)  Gn(y0) ™     h     ^       ^ 

m  =  l    „  =  I  mn  Sinh     27mnC 

(C-13) 
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APPENDIX   D 


Description  and  Explanation  of  the  IBM-650 

Computer  Program  used  to  Calculate  the 

Flux  in  the  Standard  Pile 


This  computer  program  was  written  to  calculate  the  thermal  neutron 
flux  along  the  central  axis  of  the  KSU  Standard  Pile.     The  program 
was  written  in  SOAP  II  and  floating  point.     The  object  program  and 

the  logic  diagram  are  given  in  this  section. 

r2  r2 

N     00  (n-Z+R2)     th-'  t"-1 


i  =  l    j,k=l 

z/b 


where 


2 

\    Jk/  L2  a2  b2  j         Hk 


and  erf  (t)    =    — \        e"u    du.  (D_3) 

The  program  was  set  up  to  sum  the  infinite  series    j    and    k 
by  summing  over  odd  integers    j,     holding    k    constant,   until  the 
contribution  of  a  term  of  the  series  became  less  than  a  specified 
precision.     Index    k    was  then  increased  by  two  and  the  summation 
on  index    j     repeated  until  the  precision  specification  was  met. 
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This   summation  process  was  repeated  in  such  a  manner  that  con- 
tributions from  all  possible  combinations  of  odd    j    and    k    that 
were  greater  than  a  specified  precision  were  included  in  the  sum- 
mation. 

Input  data  consisted  of  diffusion  parameters,   numerical  con- 
stants and  geometrical  specifications.     The  positions  along  the 
central  z  axis  at  which  the  flux  was  to  be  calculated  were  read  into 
drum  storage  locations  starting  with  1001.     The  program  had  a  cap- 
acity for  100  flux  positions.      The  Gaussian  ranges  were  read  into 
drum  storage  locations  0101,    0102,    and  0103.      The  corresponding 
fractional  contributions  of  each  Gaussian  term  were  read  into  drum 
locations  0104,    0105,    and  0106,    respectively.     Table  11  lists  other 
input  data  required  for  program  use. 


Table  11.     Input  data  required  for  use  of  the  IBM-650  program 
which  calculates  the  flux  in  the  standard  pile. 


Symbol 


Explanation 


Drum  storage 
location 


1  Identification  number  1  0041 

3  Identification  number  3  0969 

4  Identification  number  4  0389 
ZERO  '  0.00  0226 
ONE  1.00  '  0250 
TWO  2.00  0212 
THREE  3.00  0340 
FOUR  4.00  0407 
A  Width  of  pile  in  x  direction  (cm)  0369 
B  Width  of  pile  in  y  direction  (cm)  0469 
PI  It  0259 
SQPI  \TV  0350 
Q  Source  strength  (n/sec)  0572 
KAPP2  Reciprocal  of  diffusion  length  squared  (cm"  )  0158 
LMBDA  Transport  mean  free  path  (cm)  0025 
INDXA  Number  of  Gaussian  range  groups  0114 
INDXC  Number  of  flux  points  0202 
PREC  Specified  precision  of  last  term  in  summation  0047 


Ill 


Other  symbols  corresponding  to  terms  in  the  program  print  out 
and  the  logic  diagram  are  given  as  follows: 

FRAC        =    f . ,   the  fractional  contribution  of  the  i      Gaussian 
1 

range 

TURM       =    The  part  of  the  solution  which  depends  on  a 
particular    j    and    k 

SUM  1        =    Summation  of  TURMS 

th-i 
E  =    f.     e  (SUM  1) 


N 


i  =  l 

The  machine  yielded  an  answer  card  having  a  capacity  of  eight 
words,    a  word  being  a  ten  digit  number  and  a  sign.       Provision  was 
made  for  printing  out  partial  answers  for  checking  purposes,    thus 
making  it  necessary  to  identify  partial  answer  cards  by  the  use  of  an 
identification  number  which  was  always  printed  in  the  word  1  location 
of  the  answer  card.     Table  12  lists  the  identification  number  and  the 
answer  card  content  of  the  various  print  outs.     If  the  partial  answer 
print  outs  are  not  needed,    they  may  be  eliminated  by  use  of  the  one 
word  load  cards  listed  in  Table   13. 


Table   12.     Answer  identification  table  for  various  answer 
print  outs  for  the  standard  pile  flux. 


Word  1 

Word  2 

Word  3 

Word  4 

Word  5 

Word  6 

Word  7 

Word  8 

1 
3 
4 

J 
FRAC 
Z 

K 
SUM   1 

R 
SUM  2 

Z 

E 

TURM 

SUM   1 

F 

FLUX 

1  \l 


Table   13.      One  word  load  cards  used  lor  eliminating  unnecessary 
answer  print  outs  for   standard  pile  flux  program. 


Print  out  identification  :     Instruction  address  :        Instruction 

1  0435  60  0152  0357 

3  0263  50  0001   0183 


Table   14.     Computer  precision  study  for  standard  pile  flux. 


PREC  :        Flux  at  Convergence    Flux  at  Convergence 

z  =  5.63  cm                      time  z  a  158.08  cm  time 

10°  2259.4  n/cm2-sec  32.  361  n/cm2-sec 

10"1  2273.9                                   49  sec  32.758  17   sec 

10"2  2275.1                                    62  32.758  17 

10"3  2275.3  89  32.768  29 

10"4  2275.3  114  32.768  35 


Table  14  presents  a  study  showing  the  flux  values  as  a  function  of 
the  specified  accuracy  of  the  last  term  in  the  finite  series. 
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LDO 
LDO 

LUX 
1977 

EOO  CL 

PRINT  Z 

f   AND   FLUX 

270 

271 

272 
273 

374 

0075 
0233 
0  28  6 
0092 
OS  B  0 

81 
69 
69 

34 
69 

023  0 
03  B  6 
038  9 
197  7 
70  0  0 

0233 
0050 
009  3 
036  0 

LOO 

0453 

LDO 
STD 
LOO 
3  TO 
PCH 

1913 
LUX 
198  4 

1977 

276 
377 
378 
379 

280 

0  2  8  1 
073  5 
0  336 
0  383- 
028  7 

69 

.  ZA 
69 

04  8  2 
1.9.6  a 
033  0 

0735 

038  3 

71 

197  7 

04  77 

CONS 

1 

A 

NZC 

HLJ 

00 

0  0 

00 

0  N10 

1000 

3000 
4  000 

C  0N9 

0000 
00  00 
00  00 

NEW   Z 

1  OENT  1  F  1 

CATION 
NUMBERS 

2  S  2 
283. 
394 
-2SS 
286 

0  3  3  3 
0  33  7 
0041 
0  969. 
03B9 

a  e 

01 
00 
.00 
00 

02  11  0337 

03  8  7  03  87 
100  0  0000 

30fi-0  _Q 0-0-0 

4000   0000 

ONE 

1  0 

0000 

OOSl 

388 

0  8  5  0 

10 

00  0  0 

0051 

THREE 
FOUR 
P  1 

30 
40 
31 

0000 
00  00 
4  159 

0051 
00  51 
00  51 

390 

2.9.1 

393 

0  34  0 
.0-4  0-7- 
0259 

30 

40 
31 

00  0  0 
00  0  0 
4  15  9 

0051 
0051 
005  1 

B 
SOP  1 

17 

17 

6  330 

7  34  5 
1000 
9  3  5  D_ 
6  00  0 

00  53 
39  51 
00  57 

0  0  4  7 

00  51 

PILE  SIZE 
ROOT  OF   PI 
SOURCE   STR 

0469 

17 

623  0 

0  0  5  3 

0 

K  APP8 

16 
IS 

0573 

16 

10  0  0 

0057 

LMBD  A 

34 

MM   FR   PATH 

00  3  5 

84 

600  0 

005  1 

INOKC 

0  0 

00  00 

0014 

N  M  B  fi   OF  I 

0  203 

00 

00  0  0 

0014 

ROS 

39 

3900 

00  53 

CROUPS 

0103 

39 

390  0 

0053 

FR  AC1 

12 

BO  00 

00  5  0 

0104 

1  £ 

BO  0  0 

0050 

F  R  AC  3 

34 

1000 

00  SO 

FRACTI 0N8 

0  10  6 

34 

100  0 

005  0 

Z  1 

56 

300  0 

00  51 

Z   VALUES 

1001 

5  6 

30  0  0 

0051 

Z3 

35 

9  600 

0052 

1003 

25 

960  0 

0053 

Z5 

46 

3100 

00  58 

100  5 

46 

3100 

0052 

77 

66 

1500 

00  52 

10  0  7 

66 

1500 

0053 

29 

ZIO 

Zll 

-  2-ir» 

Z13 

a  e 

9  7 

9  40  0 
120  0 
7440 

00  52 
00  53 
00  5  3 

10  09 
1010 
1011 

86 
97 

11 

940  0 
180  0 
744  0 

0052 
0052 
0053 

IS 
—17  — 

8  080 
->40C — 

00  5  3 
Q-O-S-J 

1013 

1014 

IS 

BO  8  0 

005  3 
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APPENDIX   E 


Description  and  Explanation  of  the  IBM-650 

Computer  Program  Used  to  Calculate  the  Fast 

Point  Neutron  Source  in  a  Finite  Pile. 


This  computer  code  was  written  to  calculate  from  theoretical  con- 
siderations,  the  thermal  flux  in  a  finite  pile  due  to  a  point  source  of 
epithermal  neutrons.     The  program  was  written  in  SOAP  II  and 
floating  point.     The  object  program  and  a  logic  diagram  are  given  in 
this  section. 

The  solution  to  be  programmed  was 


CO         CO         DO 


QoF     fx)F(x0)G   (y)G   (yo)H   (z)H   (z0)e 


-(az   +02  +72)t\, 
1   m    *n       o'    th 


4.(x,y,Z)  =  )     )     y 

,    U  ,  ~   ,  D(az     +  j32   +  72     +     l/L2) 

m  =  I     n  =  l  o -1  »   m       f*n        '  0  / 


(E-l) 


where 


F     (x)     a    cos    a      x  for    m    odd  (E-2) 

m  -JST  m 


F     (x0)     =    sin  ct       x  fo 

m        *nr 


r    m    even 


and 


mir 
2a 


Fm(xo).  Gn(y)>  Gn(yo).  Ho(z0),  and  H  (z0)  are  written  in  a  similar 
manner.     The  geometry  for  this   solution  is  shown  in  Fig.    25. 

The  program  was  set  up  to  sum  the  infinite  series  by  the  process 
of  summing  over  index  o,  holding  m  and  n  constant,  until  the  con- 
tribution of  a  term  became  less  than  a  specified  precision.     Index    n 


118 


was  then  increased  by  unity  and  the  summation  on  index    o     repeated 
until  the  precision  specification  was  met.     This  summation  process 
was  repeated  in  such  a  manner  that  contributions  from  all  possible 
combinations  of  m,     n    and    o    that  were  greater  than  a  specified 
precision  were  included  in  the  summation. 

Input  data  consisted  of  various  diffusion  parameters,    numerical 
constants,    and  geometrical  specifications.     Data  were  read  into  the 
machine  with  the  program  deck  on  one-word  load  cards.     Each  one- 
word  load  card  contained  the  numerical  value  of  a  particular  constant 
and  its  specified  machine  storage  location.     Table  15  lists  the  input 
data  required  for  program  use. 


Table  15.  Input  data  required  for  use  of  the  IBM-650 
program  which  calculates  the  thermal  flux 
from  a  fast  point  source  in  a  finite  pile. 


Symbol 


Explanation 


Drum  Storage 
location 


2  Identification  number  2 

4  Identification  number  4 

6  Identification  number  6 

7  Identification  number  7 

TWO  2.00 

ONE  1 .  00 

ZERO  0.00 

PI  TT 

A  Half  width  of  pile  in  x-direction  (cm) 

B  Half  width  of  pile  in  y-direction  (cm) 

C  Half  width  of  pile  in  z-direction  (cm) 

D  Diffusion  coefficient  (cm) 

Q  Source  strength    (n/sec) 

XINTT  Initial  value  of  x    (cm) 

YINIT  Initial  value  of  y    (cm) 

ZINIT  Initial  value  of  z    (cm) 

XO  Source  location  in  x-direction  (cm) 

YO  Source  location  in  y-direction  (cm 

ZO  Source  location  in  z-direction  (cm 


0450 
0500 
0550 
0173 

0151 

0350 
0080 
0098 

0254 
0354 
0454 
0457 
0160 

0102 
0256 
0168 

0397 
0447 
0547 


Table  15.    (concl. ) 
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Symbol 


Explanation 


Drum  storage 
location 


DELX 
DELY 
DELZ 

XMAX 
YMAX 
ZMAX 

TAU 

K2 

BETSQ 
GAMSQ 

PREC 


Increment  in  x  (cm) 
Increment  in  y  (cm) 
Increment  in  z    (cm) 

Maximum  value  of  x  {cm) 
Maximum  value  of  y  (cm) 
Maximum  value  of   z    (cm) 

Fermi  age    (cm2) 

Reciprocal  of  diffusion  length  squared  (cm 

Lateral  buckling  in  y-direction    (cm       ) 
Lateral  buckling  in  z-direction    (cm       ) 

Specified  precision  of  term  in  summation 


0958 
0170 
0282 

0456 
0968 
0330 

0126 

0226 

0202 
0242 

0064 


Other  symbols  corresponding  to  terms  in  the  program  print  out 

and  the  logic  diagram  are  given  as  followsi 

FX    =    aF      (x)  F      (x0) 
m       '      m 

FY    =    b  G^  (y)  Gn  (y„) 


FZ    =    cH     (z)  H     (z0) 
o       '      o  ' 


CONSTANT  or  CONST 


-(a1     +  (32     U!) 
m      'n  o  ' 


(a*     +  02  +  72  +    k2  ) 
*    m         n         o  ' 

M,   N,   and  O    =    Summation  indices  in  the  x,   y,    and  z 

directions. 

INDXM,    INDXN,    and  INDXO    =    Summation  indices  for  test- 
ing the  oddness  or  evenness  of  M,    N,    and  O. 

The  program  calculated  the  flux  at  the  initial  values  of  x,   y,    and  z. 
The  coordinates  were  then  varied  according  to  specified  increments  in 
a  particular  manner  in  order  to  give  a  three-dimensional  mapping  of  the 
flux  in  the  region. 
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The  machine  yielded  an  answer  card  having  a  capacity  for  eight 
words,    a  word  being  a  ten  digit  number  and  a  sign.     Provision  was 
made  for  printing  out  partial  answers  for  checking  purposes,    thus 
making  it  necessary  to  identify  partial  answer  cards  by  use  of  an 
identification  number  which  was  always  printed  in  the  word  1 
location.     Table  16  lists  the  identification  number  and  the  answer  card 
content  of  the  various  print  outs.     If  the  partial  answer  print  outs  are 
not  needed  they  may  be  eliminated  by  use  of  the  one  word  load  cards 
listed  in  Table  17. 


Table  16.  Answer  identification  table  for  various  answer  print 
outs  for  thermal  flux  from  fast  point  source  in  finite 
pile. 

Word  1   :  Word  2  :  Word  3  :  Word  4  :  Word  5    :     Word  6  :  Word  7  :  Word  8 

SUM  2  SUM  3 

SUM  2  SUM  3 

SUM  2  SUM  3 

Y  Z  FLUX 


Table  17.     One-word  load  cards  used  for  eliminating  unnecessary 

answer  print  outs  for  fast  point  source  flux  in  finite  pile. 


Print  out  identification  :  Instruction  address  :  Instruction 

2  0441                         46  1077  0651 

4  0541                         46  1127  0401 

6  0591                         46  1177  0101 


2 

M 

N 

O 

SUM  1 

4 

M 

N 

O 

SUM  1 

6 

M 

N 

O 

SUM  1 

7 

XO 

YO 

ZO 

X 
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Table   18.     Computer  precision  study  for  fast  point 
source  flux  in  finite  pile. 


PREC     :     Flux  at  x  =  10.  16  cm    :     Flux  at  x  =  80.  00  cm    :  Convergence 
y  =     0.  00  y  =  80.  00  time 

z  =     0.  00  z  =  80.  00 

10°  2151.4  n/cm2-sec  2.  033  n/cm'-sec 

10"1         2153.0  2.214 

10"2         2153.1  2.207  3  min  2)    sec 

Table   18  presents  a  study  showing  the  flux  value  as  a  function  of 
the  specified  accuracy  of  the  last  term  in  the  infinite  series. 
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S 
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*       P'  ir»\>-      X 

*£  »   Y 
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-  a  o  ( 

a  >   T   * 


FLUX      FROM 


123 


BLR    1951 

8L  R    1977 
SYH   START 

SYN   X 

19  60 

198  5 
19  99 
10  00 

1 
3 

0000 
000  0 
000  0 
000  0 

00 
00 
00 
00 

00  0  0 
00  0  0 
00  0  0 
00  0  0 

0000 
0000 
00  00 
00  00 

RAM 

S  TL 

FUP 

A  A  3 

*  A  3 
A  A  3 

A  A  3 
A  A  16 

SUB 

ROUTINE 

X   FOR   C  AL 

7 

c    a 

9 
10 

00  06 
00  14 
0015 
00  32 
003  3 

20 
67 
20 
60 
39 

001  1 
0011 
0019 
001  9 
003  6 

0014 
0015 
0022 
0033 
0076 

FUP 
FAD 
FUP 

ESS 

AA  3 
AA14 
A  A  3 

SB3 

13 
1  3 

If 

000  5 
0069 
00  49 

39 
33 
39 

n 

0019 

00  7  3 

001  9 
013  3 
0019 

0069 
004  9 
0119 
009? 
0169 

FUP 

A  A  3 

A  A  11 
AA3 
A  A  10 
A  A4 

ie 

0149 

39 

0019 

0219 

FUP 
FAD 

STU 

20 
31 
32 

0199 
0269 
0349 

39 
33 
21 

0019 
037  2 
00  0  4 

026  9 
024  9 
0  00  7 

STU 
FMP 
STU 

RAU 
BM  1 

A  A  4 
A  A  4 

A  A  4 
A  A  3 
AA  5 

A  A  A  6 

34 
35 
36 
27 
36 

0054 
0057 
010  4 
0107 
006  5 

21 
39 

60 
46 

00  0  A 

00  5  7 

00  0  4 
0011 
00  16 

010  7 
0O65 
0319 

FO  V 
STU 

P.AL 

AAA 

A  A  4 

30 
3J- — 

32 

0037 
015  4 
0319 

34 
21 

65 

0004 
000  4 
000  4 

0154 

A  A  A6 

A  A  Al 

000  3 

A  A  Al  1 

34 

9998 

68  50 

34 

0233 

34   999  B 

68  50 

A  A  A13 

as 

9   1  il 

l  z  4  e 

36 

0122 

35   913  7 
17  15  6  2- 
54   302  0 

13  4  6 

A  A  A15 

54 

30  30 

0045 

3B 

□  0  39 

004  5 

EOOCR 

370 

EXIT 

COSINE 

40 

0050 

24   0053 

0056 

BM  1 

NE  C  AT 

RED  UC 

ROUTINE 

43 

0115 

46   006  8 

0  3  6  9 

BM  1 

Nf  CAT 

rnq  i  f. 

0047 

46   006  8 

0001 

R  E  D  11  C 

FSB 

au  i 

FAD 

TWOPI 

ON  EP  1 

RE  0  UC 
CO  S  1  0 

4  6 
4  7 
4B 

0369 
00  97 
0100 

33   002  1 
46   010  0 
33   020  4 

0097 
03  6  9 
0031 

RSU 

STU 

F  P  ON  £ 
TERMM 

50 
-6-1 

0039 
-  6-1-4-7- 

61   004  2 
31   000  3 

0147 
■O0-S6  — 

EO  OSR 

NEfi  A  SF 


REDUD 
-8-t  N£T— 


STO      EXIT 


FAD      T  *  OP  I 


H  E  G» 


-f  6-B — O-W-E  P  4-  ■ 

FSB      TWOPI 

FAD  ON  EP  I 
-*T-«     TM-ETA- 

RSU  800  3 

— S-T-U — T  i  BMW — 

8TU      FUNKT 


0150 
0106 
016  5 
0118 
0197 


_   0010   0013 
— a-S— 0  0  1  7 — 003  0 


010  5 
0063 
004  5 
OQ-30 
0071 


33  002  1 
_4-6  -Q-30-0 

32  020  4 
—2-4 — 0O-3-6 


—2-1 — 0002 — 0J.0  S. 


F  A  0  F  P  ON  £ 

STU  E  N  N 

RSU  TERMM 

FMP  TH  ETA 


0077 
0519 
0070 


31   0010   0063 
6  9   0  0  12  0  0  4 S 


33   00  4  3   0  4 
_2J_  "  "  " 


32  004  3 

-2-1  0  0-1-7- 

61  000  3 

39.  0-Q3-6- 

39  003  6 
—3-4 — 003  4 — 0.0-7  4 

34  00  17   0067 

31  0002   0155 


0519 
Q0.7  0 
015  7 
0.0  6  6 
013  6 


F  D  V   E  N  N 


STL   FMA  G 


0074 
_  0O.6-7- 

0155 


F  DV 

FM  A  G 

8  4 

0365 

34 

8M  1 

ENUFF 

86 

0399 

46 

00  S  3 

87 

0103 

60 

F  I1NKT 

N  F  a   H  T 

SO 

nri79 

21 

mi  n 

SIZES 

10 

0000 

0043 

8  318 

5351 

93 

10 

0000 

00  51 

31 

4159 

FP  ONE 

LOO 

I  I   Z10 

PUNCH 

96 

0156 

ROU  T  1  N  E 

97 

3_4 

98 

0030 

24 

.99 

0131 

24 

100 

1982 

8TD 

1  0  2 

00  3  4 

24 

8TD 

1984 

Z  Z  I   1 

103 

00  00 

105 

S  U  BRO  U  T 1  1 

E   106 

FAD 

BETAS 

CONSTANT 

108 

0113 

32 

FMP 

T  A  U 

" 

-10  9-  - 

111 

0  17  6 

112 

113 

0236 

20 

004  1 

00  44 

FAO   BETAS 
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FAD   CAME  116  0093  32   0046   0133 

FAD  KS  117  0123  32  0326  0253 

STU   TEMPS  118  02S3  31   0008   0061 

fit  U  T£  U.P-1- _- — JUS—-      0-0  61  60  0041  QQS-S— 

FDVTEMP2  E  X  I  N  S  130  0095  3400080303 

LOD     X  I  N  I  T  SET     II     »     2         121  1999  6901020205 

8  TO     X  A  r      l  n  i   t  l   a  L      l  -  -  0305  2410000303 

LOD      Y  I  N  I  T  VALUES  123  0303  6902560109 

STOY  124  0109  2400620365 

LDD      2IWIT 13-5 0-3-6-5 S-9 — CltS — O-t-2-1— 

S  TD      Z  C0N37  126  0121  2401240127 

SET      SUM      3         137  0127  69      0080      0133 

INOXM      A  N  U         128  0133  24      0386      0139 

U      AT      ZERO  139  0139        -     24      0092      0145 

0  N  T  3  130  0145  3400460101 

A-D-0 — Q-N-i i-3-1 0101 60     00  P  3     034  7 

TOM  132  0347  3303S00177 

133  0177  21  0092  0195 

134  0195  39  00  98  0148 

135  0148  34  0151  0201 


C0N37 

LOO 

ZERO 

8TD 

SU  M  3 

STO 

S  TD 

1  N  D  XI 

CQNT3 

■M- 

FOV 
.8-  T--U- 

FMP 
STU 
RAO 
NZE 

RAO 

4.LP-HJL 

FORM   ALPHA 

136 

020  1 

34 

02  5  4 

0  304 

C  0N11 

A  LPH2 
INDIM 
CO  NT9 

C0N11 

S  0  U  A  R  E  0 

TEST  M  FOR 
EVEN  OR   00 

138 
139 
140 

14.1 
14  2 

0111 
0  108 
0112 
0353 
0257 

39 

31 
60 
45 
60 

00  5  B 
00  5  9 
004  8 
0306 
10  0  0 

01  0  B 
0  112-  - 
0353 
03S7 
0355 

LOD 
STL 

fi  >U 
FMP 

8003 

TEMPI 
XO 

A  L  PH  » 

EOO  CR 

ALPH*   X 
FORM  COS 

14  4 
145 
146 
147 
146 

0158 
0415 
0218 
009  4. 
0251 

65 

69 
20 
60 
39 

800  3 
0218 
004  1 
03  9  7 
0058 

0415 
0050 
0094 
.03  5.1  __. 
0208 

LDD 
R1U 
FMP 
STU 
RA  U 

8002 
TEMPI 
F  X 

I  N  D  XU 

EOO  CR 

150 

046  5 

69 
6.0- 
39 

21 
60 

036  8 
8.0.0  2 
004  1 
00  9  6 
004  8 

0050 

FORM   F  X 
ADD   ONE   TO 

152 
153 
154 

0227 
0091 
034  9 

00  91 
03  4  9 
0  40  3 

C  0NT9 
C0NT8 

STU 
RA  U 
NZU 

RA  U 

Fee 

I  N  0  XM 
X  0 

C0N10 
I  N  0  XU 
ONE 

C0N12 
CON  TB 

TEST  SOR  C  E 

X  POSITION 

RESET 

156 
157 
158 
159 

160 

0277 
0306 
0213 
0318 
0453 

31 
60 

60 
33 

004  6 
D39  7 
0117 
00  4  6 
035  0 

030  1 
-0213  -- 
0318 
0453 
032  7 

C  0N10 

RA  U 
FUP 
RAL 

LDD 
STL 

ALPHA 
8  003 

TEMPI 

E008R 

FOR  M   SINE 

162 
163 
164 

165 
166 

0117 
.  0305 
035  8 
0515 
036B 

60 
39 
65 
69 
20 

1000 
005  8 
800  3 
036  8 
004  1 

0  305 
0256 
05  15 
0150-  - 
014  4 

FMP 

ALPHA 

FORM   SINE 

168 

0351 

39 
&S 

69 

005  8 
8  0.0.3 
0418 

03  0  8 

LOO 

EOO  BR 

170 

0S6  5 

0150 

FMP 

TEMPI 

Ff.au      F   » 

172 

0  377 

39 

004  1 

0141 

Rtu 

I  N  0  XM 

174 

039  9 

60 

00  4  8 

0503 

"0N1J? 

8TU 

INIKI) 

co  Nia 

FROM   INOXU 

176 

0427 

21 

00  4  8 

0301 

STO 
8  TO 
STO 

SUM2 

INDXN 

C0NT2 

1  NDXN   AND 
N   AT   ZERO 

178 
179 
180 

018  3 
0189 
024  5 

24 

24 
24 

0336   0189 
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Description  and  Explanation  of  the  IBM-650  Computer  Program 

Used  to  Calculate  the  Thermal  Flux  from  a  Thermal 

Point  Neutron  Source  in  a  Finite  Pile 


This  computer  code  was  written  to  calculate  from  theoretical  con- 
siderations,  the  thermal  flux  in  a  finite  pile  due  to  a  point  source  of 
thermal  neutrons.     The  program  was  written  in  SOAP  II  and  floating 
point.     The  object  program  and  a  logic  diagram  are  given  in  this 
section. 

The  solution  programmed  was 


*(x'y"z)=I  life     Fm<X>Fm<X°>GnWGn^ 


sinh  v      z0  sinh  y      (2c-z) 
'mn    u  'imv  ' 


v  m  m  n  n  .     ,      -, 

>mn  sinh     2y        c 


mn 


where 


(F-l) 
F     (x)  =    cos    a      x  for    m    odd  (F-2) 

m        jr  m 

sin   a      x  for    m    even 


and 


m  it 
m         2a 


F     (xn),    G   (v),    and  G   (v0)  are  written  in  a  similar  manner. 
m  a  n^'   ' 

The  geometry  for  this  solution  is  shown  in  Fig.    26. 

The  program  was  set  up  to  sum  the  infinite  series  in  the  same 
manner  as  the  infinite  series  was  summed  in  the  computer  program 
described  in  Appendix  E.        Input  data  required  for  use  of  this  pro- 
gram are  listed  in  Table   19. 
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Other  symbols  corresponding  to  terms  in  the  program  print  out 
and  the  logic  diagram  are  given  as  followsi 
FX    .    aFm(x)  Fm(x0) 
FY    =    bG   (y)  G   (y0) 


FZ 


inh  v      (2c-z)  sinh  y        z0 
'mn  '  mn 


v        sinh  2y        c 
*mn  mn 

M  and  N  =  Summation  indices  in  the  x  and  y  directions 

INDXM  and  1NDXN  =  Summation  indices  for  testing  the 
evenness  or  oddness  of  M  and  N. 


Table  19.  Input  data  required  for  use  of  the  IBM-650 
program  which  calculates  the  thermal  flux 
from  a  thermal  point  source  in  a  finite  pile. 


Symbol 


Explanation 


Drum  Storage 
location 


1  Identification  number  1 

2  Identification  number  2 

3  Identification  number  3 

4  Identification  number  4 
7  Identification  number  7 

TWO  2.00 

ONE  1 .  00 

ZERO  0.00 

PI  -n 

A  Half  width  of  pile  in  x-direction  (cm) 

B  Half  width  of  pile  in  y-direction  (cm) 

C  Half  width  of  pile  in  z-direction  (cm) 

D  Diffusion  coefficient    (cm) 

Q  Source  strength    (n/sec) 

XINIT  Initial  value  of  x    (cm) 

YINIT  Initial  value  of  y    (era) 

ZINIT  Initial  value  of  z    (cm) 

XO  Source  location  in  x-direction  (cm) 

YO  Source  location  in  y-direction  (cm) 

ZO  Source  location  in  z-direction  (cm) 


0415 
0314 
0515 
0364 
0705 

0182 
0184 
0264 
0332 

0785 
1135 
0900 
0088 
0291 

0102 
0356 
0268 

0950 
1050 
0016 
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Table   19.     (concl. ) 


Symbol 


Explanation 


Drum  Storage 
location 


XMAX  Maximum  value  of  x    (cm)  0406 

YMAX  Maximum  value  of  y    (cm)  0418 

ZMAX  Maximum  value  of  z    (cm)  0414 

DELX  Increment  in  x    (cm)  0408 

DELY  Increment  in  y    (cm)  0370 

DELZ  Increment  in  z    (cm)  0226 

K2  Reciprocal  of  diffusion  length  squared  (cm  )           0374 

TOCEE  Height  of  pile  in  z-direction    (cm)  0151 

BETSQ  Lateral  buckling  in  y-direction  (cm"  )  0376 

PREC  Specified  precision  of  last  term  in  summation         0178 


This  program  had  provision  for  printing  out  partial  answers  for 
checking  purposes.     Table  20  lists  the  identification  number  and  the 
answer  card  content  of  the  various  print  outs.     If  the  partial  answer 
print  outs  are  not  needed  they  may  be  eliminated  by  use  of  the  one 
word  load  cards  listed  in  Table  21.     Table  22  presents  a  study  show- 
ing the  flux  value  as  a  function  of  the  specified  accuracy  of  the  last 
term  in  the  finite  series. 


Table  20.     Answer  identification  table  for  various  answer 
print  outs  for  thermal  flux  from  a  thermal 
point  source  in  a  finite  pile. 


Word  I 

:  Word  2  : 

Word 

3  : 

Word 

4  : 

Word  5 

:  Word  6  : 

Word  7  : 

Word  8 

1 

M 

N 

FY 

FZ 

SUM  1 

SUM  2 

2 

M 

N 

FY 

FZ 

SUM  1 

SUM  2 

3 

M 

N 

FX 

FZ 

SUM  1 

SUM  2 

4 

M 

N 

FX 

FZ 

SUM  1 

SUM  2 

7 

XO 

YO 

ZO 

X 

Y 

Z 

FLUX 

Table   21.     One  word  load  cards  used  for  eliminating 
unnecessary  answer  print  outs  for  thermal 
point  source  in  a  finite  pile. 
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Print  out  identification  :     Instruction  Address  : 


Instruction 


0605 
0258 
0655 
0308 


46  0258  1059 
69  0184  0587 
46  0308  0735 
60  0220  0425 


Table  22.       Computer  precision  study  for  thermal  point 
source  flux  in  finite  pile. 


PREC        Flux  at  x  =  0 
y  =  0 
and  z  =  20.  32  cm 


Convergence      Flux  at  x  =  0  Convergence 

time  y  =  0  time 

and  z  =  81.72  cm 


10°  4449.  1  n/cm2-sec 

10"1  5202.2 

10"2  5339.7 

10"3  5363.3 


approx.    5 
min. 


107.  31  n/cm2-sec 

126. 16 

129.07 

129.47 


approx. 
30  sec. 
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START 


CON  26 

/Set    SUM  2\ 
INDXM     a 
\   ot    ZERO 


Read     in        \ 
'rogram    and 
Constants 


y-* 


CONT    5 


Subtract    ONE 
from     INDXM 


CONT    7 


Add    ONE    to 

N.    Form     p„ 

a    H 


Subtroct   ONE 
from    INOXM 


Calculote 

FV-cos/9„>cosfl„»0 

Add    ONE    to 

INDXN 


CON    10 


Calculate 
FY.sin/3ny  sin/S„y0 
Subtract   ONE 
from     INOXN 


CON    12 


Calculate 


To  CONT  7 


To  CONT    I 


To  CON   26 


Set  X  8  XINITN 
-   XINIT  + 

DELX.     Set    Y 
Y    +     DELY 


To  CON    26 


LOGIC   DIAGRAM- APPENDIX  F 
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33 

040  8 

17  8  5 

F  SB 

X  MA  X 

471 

0553 

33 

040  6 

06  63 

C  0  r.  ;■  ■?.              R«U 

X  1  N  1  T 

RESET   X 

473 

068  7 

6  0 

010  2 

0857 

STU 
ATM 

X   1  N  1  T 

475 

1835 

B] 

010  3  0605 

RAU 

CHANGE   V 

477 

060  3 

60 

0313 

0517 

STU 

Y 

MAXIMUM   Y 

479 

1097 

81 

0313 

0615 

BH  1 

CO  N36 

C0N2S 

481 

044  5 

46 

036  1 

0649 

06  9  9-    Q&9  9— 
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1 

00 

1000 

0000 

1 OE  NT  1  F  1 

483 

04  15 

00 

100  0 

000  0 

3 

00 
00 

00 

3000 
4000 
7000 

0000 

00  00 
0000 

NUMBER  8 

485 

0515 

00 

3000  0000 

T 

487 

0705 

00 

7000 

0000 

ONE 

10 

0000 

0051 

489 

0184 

10 

000  0 

0051 

P  1 

31 

4161 

0051 

491 

0333 

31 

416  1 

0051 

B 

88 

1100 

0053 

PILE   SIZE 

1135 

88 

110  0 

0052 

a 

82 

0000 

0050 

DIFF   COEFF 

0088 

82 

000  0 

0050 

X  I  N  1  T 

00 

0000 

0000 

INITIAL  X 

0102 

00 

oo  o  a 

0000 

Z  t  N  1  T 

96 

2700 

0053 

INITIAL  Z 

0268 

98 

270  t 

0052 

YO 

00 

0000 

0000 

source  pos 

1050 

00 

0000 

0000 

KM  AX 

se 

1  100 

0053 

MAXIMUM 

0406 

68 

1100 

0052 

OELY 

20 

3300 

0052 

X   Y   AND   Z 

0  370 

20 

330  0 

005  2 

OELZ 

10 

1600 

00  52 

02  26 

10 

160  0 

0053 

K3 

35 

9250 

0047 

KAPPA   SOO 

0374 

35 

925  6 

0047 

BE  T80 

7840 

0047 

Y   BUCKL 1 NG 

0376 

31 

764  0 

0047 

PRE  0  1  8  1  ON 

0176 — 
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APPENDIX   G 


Description  and  Explanation  of  the  IBM-650 
Computer  Program  Used  for  Fitting  Empir- 
ically Experimental  Data  with  the  Sum  of 
Several  Terms  of  Gaussian  Form 


The  computer  code  was  written  to  fit  empirically  the  experi- 
mental indium  resonance  activity  with  the  sum  of  several  Gaussian 
terms  having  adjustable  parameters.  The  program  was  written  in 
SOAP  II  and  floating  point  form.  The  object  program  is  listed  and 
the  logic  diagram  is  given  in  this  section. 

This  method  of  fit  was  not  exact.     There  was  no  guarantee  that 
the  result  of  this  fit  would  be  the  best;    but  in  a  practical  approach 
to  the  problem  of  fitting  a  non-linear  function,   the  result  appeared 
to  be  adequate. 

The  criterion  that  the  machine  inspected  was  that  the  sum  of  the 
weighted  squares  of  the  residuals  between  the  experimental  data  and 
calculated  values  should  be  made  as  small  as  possible.     This  was  done 
through  a  trial  and  error  process.     Each  of  the  fitting  parameters 
was  varied  in  turn  by  a  specified  increment,  holding  all  other  para- 
meters constant,   until  such  time  that  a  specified  increment  could 
make  no  further  reduction  in  the  sum  of  the  weighted  squares  of  the 
residuals.     This  parameter  was  then  stored  as  the  best  available  esti- 
mate of  the  particular  empirical  parameter.     When  none  of  the  para- 
meters could  be  varied  by  the  specified  increment  to  give  a  smaller 
sum  of  the  weighted  squares  of  the  residuals,   the  increments  were 
refined  and  the  trial  and  error  process  was  repeated  with  the  refined 
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increments.     This  procedure  continued  until  the  increments  were  less 
than  a  specified  precision. 

The  second  portion  of  the  program  calculated  the  fractional  con- 
tributions of  each  range  group  according  to  Eq.    16  of  the  theory. 

The  data  were  empirically  fit  with  a  function  of  the  form  of  Eq. 
14  of  the  theory 

N 


FLUX 


.    =    )    B.  e      J  ,      N  g  10  (G-l) 


i  =  l 

As  stated  above,    the  best  fit  criterion  was  that 


M  /  N  x2 


ERROR     a    )      — !- I     )      B.   e        ■>      j     -  FLUX.]       ,    MS  20  (G-2) 

be  a  minimum. 

The  program  required,    in  addition  to  the  experimental  data  and 
their  respective  positions,   initial  estimates  for  a  specific  number 
of  parameters.     The  program  had  a  capacity  for  up  to  20  data  points 
and  10  terms  in  the  summation  of  Eq.    G-l.     These  input  data  were 
read  into  the  machine  along  with  the  program  deck  on  one -word  load 
cards.     Each  one -word  load  card  contained  a  particular  constant  or 
initial  value  and  its  specified  storage  location.     Table  23  lists  the 
various  input  data  needed  for  this  program. 


Table  23.     Input  data  required  for  use  of  the  IBM-650 

program  which  empirically  fits  experimental 
data  with  several  terms  of  Gaussian  form. 
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Symbol 


Explanation 


Storage 
location 


ZERO 
FPONE 
FPTWO 
FOUR 

TEN 

ONE 

EIGHT 

ONEHD 

PI 

SQPI 

PREC 

Al 
A2 

DELB 
DELR 

INDXA 
INDXB 

Z 

FLUX 
W 


0.  00 
1.00 
2.  00 
4.  00 
10.00 

Index  number  1 
Index  number  8 
Desired  precision  of  parameter  increments 


Desired  precision  in  term  of  series  summation 

Pile  width  in  x-direction    (cm) 
Pile  width  in  y-direction    (cm) 

Initial  increment  of  the  amplitude,     B. 

Initial  increment  of  the  Gaussian  range,    r.(cm) 

Number  of  data  points 

Number  of  Gaussian  terms  in  empirical  fit 

Position  of  j       data  point    (cm) 
Experimental  data  at  j  ,    point 
Weighting  function  of  j       piece  of  data 
Initial  estimate  of  B. 


Initial  estimate  of  r.     (cm) 


0059 
0096 
0260 
0743 
0436 

0292 
0124 
0354 
0290 
0164 
0126 

0405 
0211 

0028 
0198 

0077 
0162 

12(00+j) 
12(20+j) 
17(50+j) 
13(00+i) 
13(20+i) 


The  machine  yielded  an  answer  card  having  a  capacity  of  8  words, 
a  word  being  a  ten  digit  number  and  a  sign.     For  the  first  answer 
punching  operation  the  machine  printed  out  the  initial  estimates  of  the 
fitting  parameters  on  as  many  cards  as  was  necessary  to  accomodate 
them.     The  amplitude  and  range  of  the  first  term  were  printed  out  in 
word  locations  1  and  2,   respectively.     The  amplitude  and  range  of  the 
second  term  were  printed  in  word  locations   3  and  4,    respectively. 
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After  a  card  was  filled  to  its  8  word  capacity  it  was  punched  and  a 
new  card  began  to  be  filled.     This  procedure  continued  until  all 
fitting  parameters  had  been  punched  out.     Then  a  separate  card  was 
punched  giving  in  the  word  8  location,   the  value  of  the  sum  of  the 
weighted  squares  of  the  residuals  between  the  experimental  data  and 
the  calculated  value.     Subsequent  improvements  in  the  parameters 
were  printed  out  according  to  this  procedure  after  each  cycle  of  try- 
ing to  vary  each  parameter.     When  the  fitting  parameters  could  not 
be  further  improved  with  the  most  refined  increment  specified,    the 
punching  out  of  the  best  fit  parameters  and  their  sum  of  the  weighted 
squares  of  the  residuals  took  place  according  to  the  procedure  pre- 
viously described.     At  this  time  the  residual  at  each  point  was  punched 
out  along  with  the  position  of  the  point.     The  position  of  the  last  point 
and  the  associated  residual  were  loaded  into  word  locations  1  and  2, 
respectively.     The  position  and  the  residual  of  the  next  to  last  point 
were  loaded  into  word  locations  3  and  4,    respectively.     This  pro- 
cedure continued  until  all  the  residuals  and  their  identifying  locations 
had  been  punched  on  as  many  cards  as  was  necessary. 

The  machine  then  calculated  the  fractional  contributions  of  each 
range  group  and  printed  it  out  in  the  word  5  location  along  with  the 
respective  Gaussian  range  in  the  word  7  location. 
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SYN 
SYN 
8YN 
SYN  _ 

SYN-  #  ye  st 


POSITION 


OOOO 
000  0 
0000 
0000 
0000 
000  0 
0000 
0000 

flo-a-o— 


000  0 
00  0  0 
000  0 
00  0  0 
000  0 
000  0 
000ft- 


0000 
0000 

oo-o-o— 


R 

1280 
1390 
1300 
1310 
1330 

RANGE 

10 

0000 

00 

000  0 

0000 

a  i  N  i  T 

INITIAL 

12 

000  0 

00 

000  0 

000  0 

TEMP 

14 

000  0 

00 

000  0 

0000 

C 

F  R  A  C 

0  U  Ot 

A  A  A  1 

17  30 

CONSTANT 

1  6 

0000 

0  0 

00  0  0 
000  0 
00  0  0 

0000 

1740 

1  B 

0000 

00 

0000 

E   TO   X 

30 

000  0 

84 

000  3 

00  06 

R  AU 

A  A  A  2 

1      FOR  C.   Al  C    23 

00  14 

67 
20 

0011 
Qfll  9 

0015 

no?,  p. 

RA  U 

A  A  A  3 

24 
35 
36 

0022 
0  02  3 

007  6 

39 
32 

002  6 

0076 

FAD 

A  A  A  15 

003  9 
001  9 

OOOS 

0  fl  fi  <j 

FAD 

A  A  A  14 

36 

0069 
0049 

32   00  7  2 

30  nn  1  o 

004  9 

m  1  9 

FAD 

A  AA13 

30 

0119 
0099 

32 

39 

013  3 
00  19 

009  9 
01  «9 

FAO 

A  A  A13 

32 

0169 

32 

017  3 

014  9 

FAD 

A  A  A  11 

34 

0  319 

32 

03  3  2 

019  9 

FAD 

A  A  A10 

36 

0369 

33 

037  3 

024  9 

FMP 

A  A  A  4 

36 

0007 

39 

00  0  4 

0054 

FMP 

A  A  A  4 

4  0 

0057 

39 

00  0  4 

0104 

RA  U 

A  A  A  3 

42 

010  7 

60 

00  11 

006  5 

AA  AS 

RIU 

A  A  A  10 
1  At  4 

44 

00  16 

60 

037  3 

nn  n  4 

002  7 

1.4 

STU 

A  A  A  4 

A  A  A6 

46 

0154 

21 

00  0  4 

0  3  19 

A  A  A10 

10 

0000 

0051 

4  S 

0  372 

10 

000  0 

005  1 

AA  A13 
A  A  A13 
A  A  A14 

31 
35 
17 

3575 
9  137 
1562 

83*9 
1248 
0047 

50 
52 

0172 
012  3 
0073 

31 
35 
17 

357  5 
913  7 

156  3 

834  9 
124  8 
004  7 

A  A  A16 

69    0  600 

004  4 

54 

0026 

69 

060  0 

004  4 

LOO 

Z  Z  Z10 

PUNCH   SAND 

56 

00  56 

69 

00  0  9 

0013 

6T0 

1976 

58 

0030 

24 

197  8 

00  3  1 

STO 

I960 

60 

0033 

34 

198  0 

0033 

STO 

1983 

63 

00  34 

34 

198  3 

003  5 

fl 

8  TO 

1984 

Z2Z1 

6  4 

00  36 

198  4 

00  5  3 

TEMPP 

3TD 

EG  8  1  T 

ROUT  1  NE 

66 

0100 

24 

010  3 

010  6 

FMP 
FOJL 
FDV 

Z 

J 

RANGE 

68 

00  55 

39 

320  0 

015  0 

R 

e 

PART   OF 

70 

0060 

34 

528  0 

0130 

LOO 

EOOE  A 

73 
73 

0037 

69 

_2fl_ 

004  0  0000 

TURMU 

BTO 

EX  X  1  T 

RO  UT  t  N  E 
T  n   r  r>  □  u 

74 
7S 

0200 
0116 

24 

0153 

7.3  3  n 

0156 

I 

1 

TERM 

OF 

76 

0025 

39 

536  0 

0010 

EXIT 

78 

0250 

030  3 

020  6 

RA  C 

aooi 

TO 

FORM 

80 

0062 

88 

800  1 

0  068 

LOO 

1  N  D  XA 

HE  1 

GHTEO 

82 

00  34 

69 

007  7 

018  0 

R  E  AOY 

LDO 

1  N  0  XC 

SOU 

ARE  0 

64 

0086 

69 

00  3  9 

004  2 

LDO 

ZERO 

nnnuE 

86 

004  8 

69 

00  5  9 

0112 

DOONE 

R  A  U 

TERM 

C 

SB 

one 

60 

752  0 

007  5 

STO 

A  CCUM 

90 

0041 

21 

0115 

0  1  (',  P 

SXC 

B001 

92 

0230 

59 

800  1 

0136 

R  EPET 

BMC 

D  1  0  1  T 

DOONE 

94 

0089 

j  i 

00  9  0 

0118 

FSB 

FLUX 

I 

96 

0  369 

33 

323  0 

0  0  A  ■? 

FMP 

G 

A 

98 

004  3 

39 

324  0 

014  0 

EXIT 

ami*  in 

104 

0186 

40 

00  8  6 

020  3 

LDO 

EOO  CL 

To   PR  1 

T 

106 

0256 

69 

010  9 

005  0 

RA  B 

8  001 

REB  1  T 

108 

016  5 

62 

800  1 

0071 

R3A 

8  001 

FILL 

110 

0127 

81   80  0  1 

0063 

STO 

196  S 

A 

112 

00  7  3 

34   39B  5 

003  B 

BTO 

1986 

A 

114 

0093 

2  4   3'9  8  6 

013  9 

NEXTB 

AXA 
N  I  B 
NZ  A 

0003 
NEXTB 
FILL 

1  HO  P  E 

C  A  R  DZ 

116 
117 

lie 

004  5 
0001 
0204 

50   00  0  2 
42   030  4 
40   00  8  3 

000  1 
0105 
000  fl 

PCH 
LOD 

—tee 

S  TO 


EOOCL 
EOOCL 


119 
120 
131 


010  5 

033  7 

-S-3-3-0— 

0336 


007  1  ( 
197  7  ( 
03  3  0  ( 
■01-3  3-  -i 


PCM 
3T0 
LOD 
S  TD 

1977 
T  1  XE 

B 
BREST 

B 
B 

[XET 

ROUT  1  N 
TO  SET 
BEST 

008  7 
0  400 
0  30  6 
00  13 

71 
34 
69 

24 

197  7 
030  3 

536  0 

537  0 

PRIME 

E 

136 
137 

12B 

0306 
0013 
01  33 

-START 

to aoy — 

8T0 
LOO 
STD 

RBE8T 
ERROR 
ERHIN 

B 

TIKE 

PARAMETERS 

130 
131 

133 

0  1  B  3 
014  3 
0  17  4 

24 
69 

24 

S3  9  0 
00  3  1 
013  3 

0  14  3 
017  4 
0  30  3 

HUB 
— UH>- 
8TD 
STD 
LDD 
STD 
STO 

B  001 

LOADY 

134 

0  2  15 

63 

BOOl 

012  1 

B 

BBE3T 

R  1  N  1  T 

R 

P.  B  E  ST 

B 
B 
B 
B 
B 

PAR  AM  E 

at  mi 

E  S  T  )  II  A 

TER  S 

136 

0  353 

34 
34 
69 

24 
24 

53  6  0 
53  7  0 
5310 
538  0 
529  0 

0063 

TE  S 

138 
13  9- 
140 

0173 
0113 
0233 

0113 

0  19  3 

-SETUP 

NZB 

Rill 
R  *  » 
LDD 
R*  B 

LO  A  DY 

SETUP 

148 

0299 

43 

012  1 

040  3 

BOOl 
INOXB 
BOOl 

144 

0081 

60 

BOOl 

0137 

DETERM 

N  E 

146 

82 

BOOl 

017  1 

C  AL  C  A 

TEMfP 
T  UR  MM 


148 

004  4 

30 
88 
69 

00  3  9 
BOOl 
00  5  1 

0093 

150 

0098 

010  0 

1S2 

025  4 

59 

000  1 

006  0 

NZ  A 

CAL  CA 

C  AL  C  B 

COMB  1  N 

154 

0016 

oo  g  a 

002  0 

R  AA 

BOOl 

0001 

CA  LCA 

FILUP 

FBO 

OF   R   AND   B 

156 
157 
158 

0380 
0386 
014  2 

BO 
53 
42 

aoo  i 

00  0  1 
00  9  B 

03  B  6 

NZB 

00  4  6 

LDD 

—BA-B- 

LD  0 


161 

162 

.  lfi.3  - 

164 


0002 
0155 
0315 
0221 


69      0003   0400 
6  9   015  5   0  3S.Q- 
0315 
022  1.. 


69 


016  2 
80  0  1 

00  3  9 


ZERO 
i  NDXF 
BBEST 


Bfi   ROOT 01  4 


00  5  9 
34  036  5 
60  537  0 
3  3  00  2  B. 

526  0 


i.  n  n 

S  TO 

CT  EUP 

172 

0  4  19 

2  4 

032  3 

0175 

8T0 

6TEMP 

174 

0131 

24 

00  8  4 

0187 

R  A  A 

8001 

S  U  BB2 

i  ilB.il     UL* 

176 

0  430 

BO 

BOOl 

0336 

SXC 

0001 

TERM 

178 

01S9 

59 

00  0  1 

009  5 

HZ  A 

8  U  BBS 

5  11  BS  3 

OET  ER  M  1  NE 

180 

0101 

40 

0  3  3  6 

0355 

RAU 

FSB 
BM  1 

Efl  M  1  N 

ERROR 
B  A  D  OB 

GOO  DB 

TEST  ERROR 

182 
183 
184 

0058 
0337 
014  7 

60 
33 
46 

013  3 

003  1 

04  5  0 

0337 
0147 
0151 

R  AB 

8001 

t  MPR  0  V  E  0   B 

186 

02B7 

BS 

80  0  1 

024  3 

RAC 

8  001 

181 

0235 

88 

BOOl 

0181 

F  SB 

BM  t 

B         B 
BP  LUS 

BU  1  M  3 

190 

0375 

3  ) 

536  0 

0337 

BP  LU  B 

LDO 

PRIME 

BET  B 

192 

0  190 

69 

039  3 

040  0 

STD 

1  NOXF 

194 

0399 

84 

0  36  5 

036  B 

R  A  U 

BBEST  8 

INCREASE  B 

196 

0  23  4 

60   527  0 

0  335 

B  M  1    Hi, 

LOO 

PRIME 

SET  B 

198 

009  1 

69 

00  9  4 

04  00 

STD 

INDKf 

SOO 

0262 

24 

036  5 

0318 

STD 

1  N  0  XO 

30  3 

0  449 

14 

03  7  1 

0374 

FSB 

OELB 

SUBB1 

204 

0  37  5 

33 

002  6 

0  30  5 

R  AB 

BOOl 

IMPROVED   B 

206 

0  387 

a  a 

BOOl 

0343 

NZ  U  LOOPR 

LETUP 

209 

0469 

44 

0323   0334 

R  AB 

8001 

210 

0  437 

es 

800  1 

0393 

FSB 

BU  t 

B       B 

M.O.fiE-3 

i  rtnon 

313 

51  T 

0425 

niiflT 

33 

5  26  0 

04  8  7 

LDD 
fi  AB 
LOD 
RAC 

RAU 
FSR 


8T£M 

B  001 

CTEMP 

BOOl 

BBEST 

nF  L  R 

ZERO 
8  TD   I  N  OX.F 
LDO   BBEST 
8TD  B 
RAU   RBEST 

_F_i-D "FIB 

STU     R 
LDO     C  T  E  MP 
RAC         8001 
LD  0      8  T  EUP 
R A  B  S001 


HZ  A aitifna 

RAU      RBEST 


215- 

216 

317 


DECREASE      8      21 


I  N  C  R  F.  i  .S-E B 325 


69  00  8  4 

83  BOOl 

69  032  2 

_B_8_  8  0  0  1 

60  527  0 


32 QJULa 0  5  7  5 

21      5280      0283 
69      032  3      Q62_5_ 
SB      aOOl      0281 
6-9.      008  4      OS  8-3 
83      BOOl      0493 


8001 

023  9 

024  2 
00  01 
00  0  1 


0386 
0  1  0-0- 
020  0 
0  2  48— 
0304 


143 


144 


ERROR 

241 
242 

0061 
0064 

69 
60 

006  4 

013  3 

0250 

RAU 

ER  M  1  N 

.Q19-7 

GOODS 

BM  1 
LOO 
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APPENDIX  H 


Consideration  of  Perturbations  and  Corrections  Encountered 
in  Experiments  Using  Finite  Foils  and  Sources 


In  order  to  determine  indium  foil  activations  free  from  self- 
perturbations  or  experimental  limitations,    it  is  possible  to  consider 
the  following: 

(a)  activation  of  indium  by  higher  energy  resonance  neutrons 

(b)  absorption  of  epithermal  neutrons  in  cadmium  filter 

(c)  finite  size  of  source  and  detector 

(d)  difference  in  foil  weights  or  size 

(e)  unequal  activation  of  two  sides  of  foil 

(f)  flux  depression  in  vicinity  of  foil 

(g)  self-shielding  of  the  center  of  the  foil  by  its  outer  layers 
(h)    mutual  shadowing  of  neighboring  foils 

(i)    hardening  of  the  neutron  spectrum  through  the  foil 
(j)    obliquity  of  foil  with  respect  to  the  source. 
If  there  is  need  for  an  estimate  of  the  absolute  flux  in  medium  or 
physical  situation  that  differs  from  the  standardizing  medium  or  the 
standardizing  physical  situation,   these  factors  should  be  considered. 
Factors  (a)  through  (e)  have  been  considered  in  this  pile  standardi- 
zation. 
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Activation  of  Indium  by  Higher  Energy  Resonance  Neutrons 

This  section  considers  the  correction  factor  that  must  be  applied 
to  indium  foil  measurements  made  near  high  energy  neutron  sources, 
because  of  activation  of  the  indium  by  resonance  neutrons  about  1.44 
ev.    (19,    41). 

According  to  Hill,    Roberts  and  McCammon,    (19),    at  points  close 
to  a  source  of  high  energy  neutrons  in  graphite,    a  good  part  of  the 
foil  activation  is  due  to  absorption  of  neutrons  above  an  energy  of 
1.44  ev.     Farther  from  the  source,    the  perturbation  due  to  higher 
resonances  becomes  less  so  that  this  correction  factor  becomes 
negligible. 

The  experimental  technique  for  determining  the  correction  factor 
to  be  applied  to  the  measured  activity  of  the  cadmium  covered  indium 
foils  consists  of  comparing  the  activity  of  indium  foils  covered  by 
both  cadmium  and  indium  with  the  activity  of  foils  covered  by  cadmium. 
Since  the  resonances  of  indium  115  above  1 .  44  ev  are  weak  compared 
with  the  resonance  at  1.44  ev,   the  self- shielding  of  the  high  energy 
resonances  is  small  compared  with  that  of  the   1 .  44  ev  resonance. 
Consequently,    covering  the  foils  with  indium  will  have  little  effect  on 
the  activation  by  neutrons  of  higher  energy,   but  will  decrease  the 
activation  at  1.44  ev  by  a  large  factor. 

At  distances  less  than    \Tt      from  the  source  (41),    the  total  acti- 
vation of  cadmium  covered  foils  is 

CdlnCd  1.44  ev       "higher  energy.  (H-l) 
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At  distances  of  two  or  three  times    \[r~  from  the  source,   the 
activation  by  high  energy  neutrons  is  negligible  so  that  the  activa- 
tion is 

A_,T    L,A    x    A.     ..  (H-Z) 

CdlnCd  1 .  44  ev  s         ' 

For  the  case  in  which  the  foils  are  covered  with  both  indium 
and  cadmium,   Eqs.   (H-l)  and  (H-Z)  can  be  written,    respectively 

ACdInInInCd     =    go    A  1 .  44  ev    +    A higher  energy  ^H"3^ 

ACdInInInCd    =    go  A  1 .  44  ev  (H~4) 

where    g0    is  a  shielding  factor  for  the  indium  foils  at  the 
1.44  ev  resonance.     Equations  (H-Z)  and  (H-4)  maybe  used 
to  obtain    g0    and  then  Eqs.    (H-l)  and  (H-3)  maybe  solved  simultan- 
eously to  obtain     A,  .    ,  from  known  measured  values  of 
higher  energy 

ACdInCd    and    ACdInCd- 

Additional  experimental  measurements  were  made  by  enclos- 
ing the  detecting  indium  foil  on  each  side  with  similar  indium 
foils  in  addition  to  the  usual  cadmium  cover.     Cadmium  covered 
indium  foil  measurements  made  by  Steichen  (39)  were  used  for  the 
unshielded  measurements  after  verifying  that  Steichen's  measure- 
ments could  be  reproduced. 

Table  Z4  lists  the  experimental  results  of  foil  measurements 
and  the  ratio  of  the  indium  foil  activities  when  covered  with  indium 
plus  cadmium  to  the  activities  when  covered  only  with  cadmium. 
Fig.    Z7  is  a  plot  of  this  ratio. 
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By  assuming  that  the  ratio  will  approach  the  constant  shielding 
factor,     g0,     at  points  away  from  the  source,    a  line  representing 
this  shielding  factor  has  been  approximated.     With  this  approxima- 
tion,   a  correction  factor  has  been  determined  and  plotted  in 
Fig.    28  along  with  the  results  of  other  work  done  with  this  cor- 
rection factor  (19,41). 

Absorption  of  Epithermal  Neutrons  in  the  Cadmium  Filter 

This  correction  factor,   designated  by    F     ,,     considers  the 

■  '        cd 

absorption  of  epithermal  neutrons  in  the  cadmium  foil  covering. 
It  is  a  function  of  both  the  thickness  of  the  cadmium  covering  and 
the  thickness  of  the  foil.     For  the  particular  foils  and  covers  used 
in  this  work,     F         was   1.16.     The  graph  from  which  it  was  ob- 
tained is  available  in  many  references,    (2,    18,    28,    40). 

This  factor  was  measured  by  Tittle  (40)  by  a  technique  that 
used  successively  different  thicknesses  of  cadmium  covering  for 
a  particular  foil  thickness,    the  correction  factor  being  found  from 
an  extrapolation  to  zero  cadmium  thickness. 

Finite  Size  of  Source  and  Detector 

To  correct  for  the  finite  sizes  of  the  source  and  detector,   the 
following  equation  may  be  used  for  a  circular  source  of  radius  a, 
and  a  circular  detector  of  radius  a2, (41). 
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A(ro).Am(ro)     .     .(-£^L,C*A)  (H-5) 

r=o 

A     (r0)    is  the  observed  activity  at  the  distance    r0    between 
the  centers  of  the  source  and  detector.     A(r0)  is  the  ideal  activity 
of  a  point  source. 

Steichen  (39)  has  determined  this  correction  factor  for  resonance 
neutron  measurements.     It  is  given  in  Fig.    29.     It  was  applied  to  the 
resonance  measurements  made  in  this  work.     This  factor  would  vary 
with  the  diffusing  medium. 

Variation  in  Foil  Weights  or  Sizes 

Graves  and  Froman  (16)  point  out  that  for  thick  foils,   variation 
in  weight  has  a  small  effect  so  long  as  the  foil  area  is  constant  from 
foil  to  foil.     They  state  that  only  the  activity  of  a  thin  foil  is  pro- 
portional to  the  mass  of  the  foil,   and  that  a  correction  for  variation 
of  foil  mass  is  difficult  to  make  for  a  foil  of  sufficient  thickness  to 
be  usable.    They  feel  that  it  is  better  to  control  the  manufacture  of 
the  foils  such  that  the  variation  in  weight  from  foil  to  foil  is  very  small. 

In  the  fitting  of  the  resonance  activations  described  in  the  analyt- 
ical portion  of  the  work,   it  was  found  that  consideration  of  the  varia- 
tion of  the  foil  weights  improved  the  results  of  the  resonance  data  fit 
considerably.     The  average  foil  weight  was  0.4592  grams  with  a  maxi- 
mum deviation  of  one  per  cent  by  any  foil. 
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Fig.  29.     Source  and  detector  geometry   correction  factors 
for   the  single    source  and  four   source  resonance 
distributions. 
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Unequal  Activation  of  Two  Sides  of  Foil 

Since  foils  are  seldom  thin  for  neutrons  or  for  the  electrons 
produced  by  induced  activation,   the  activity  measured  on  different 
sides  of  the  foil  may  be  different.     This  difference  depends  upon 
the  flux  gradient  in  the  neighborhood  of  the  foil,     (16).       It  is 
generally  agreed  (2,    16)  that  whenever  the  diffusion  approxima- 
tion is  adequate  to  describe  the  neutron  flux,   the  average  of  the 
measurements  on  both  sides  of  a  foil  is  proportional  to  the  flux. 

All  foil  data  reported  in  this  work  is  the  average  of  measure- 
ments made  on  both  sides  of  the  foils. 

Flux  Depression  in  the  Vicinity  of  the  Foil 

A  flux  detecting  foil  absorbs  neutrons  which  would  otherwise  be 
able  to  scatter  back  and  bolster  the  flux  at  the  foil  position.     The 
magnitude  of  this  flux  depression  depends  upon  the  properties  of  the 
foil,   the  angular  distribution  of  the  neutrons  in  the  foil  vicinity, 
and  the  mean  free  path  of  neutrons  in  the  surrounding  medium  (30). 

This  factor  need  not  be  considered  when  measurements  with  the 
standardized  detector  are  made  in  the  standardizing  medium.     If  the 
standardized  detector  is  used  in  other  mediums,   this  factor  may  be 
easily  determined  from  formulas  given  in  various  references, 
(28,    31,    34). 

This  correction  for  flux  depression  need  not  be  considered  for 
measurements  of  slowing  down  density  at  indium  resonance  energy. 
Absorption  of  an  epithermal  neutron  does  not  remove  a  neutron  which 
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could  be  scattered  back  to  bolster  the  fluxi   since  in  order  for  the 
neutron  to  be  scattered  back  to  the  detector  position,   the  neutron 
must  suffer  an  energy  degrading  collision.     Thus  the  neutron  may 
not  have  the  resonance  energy  and  would  not  be  eligible  to  bolster 
the  resonance  energy  flux. 

Self-protection  of  the  Foil  by  Its  Outer  Layers 

'  This  consideration  takes  into  account  the  fact  that  atoms  in  in- 
ternal layers  of  a  thick  detector  are  exposed  to  a  neutron  flux  smaller 
than  that  impinging  on  the  surface  of  the  foil.     The  self -protection 
ratio  is  the  ratio  of  the  gross  activation  of  the  foil  to  the  quantity  cal- 
culated by  multiplying  the  macroscopic  cross  section  for  activation 
by  the  surface  flux  (31).     Determination  of  this  self-protection  ratio 
is  outlined  in  reference  (31). 

This  factor  does  not  depend  upon  foil  position  and  may  be  con- 
sidered a  constant  for  any  medium  and  position.     It  is  a  function  of 
only  the  foil  material  and  thickness,    and  the  nature  of  the  flux.     It 
would  be  necessary  to  apply  this  correction  factor  if  a  collimated  flux 
were  being  measured.     For  the  assumed  isotropic  flux  in  this  work, 
this  correction  factor  has  not  been  considered. 

Mutual  Shadowing  of  Neighboring  Foils 

It  is  possible  to  irradiate  resonance  detectors  near  enough  to  each 
other  to  cause  an  effect  by  one  foil  upon  the  other.     Approximate  values 
for  the  fractional  loss  of  resonance  activations  due  to  the  presence  of 
a  neighboring  foil  are  given  in  reference  (30).     For  the  particular 
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foils  and  foil  spacings  used  in  this  work,    the  effect  of  one  foil  on  the 
other  was  estimated  to  be  negligible. 

Hardening  of  the  Neutron  Spectrum  through  the  Foil 

The  correction  for  the  hardening  of  the  spectrum  arises  from  the 
fact  that  neutrons  of  lower  energy  are  preferentially  absorbed  in  the 
first  few  layers  of  a  thick  foil.     Thus  the  neutrons  reaching  the  inter- 
nal layers  of  the  detector  will  be  of  higher  average  energy  and  will 
be  subject  to  different  capture  cross  sections. 

Amaldi  (2)  and  Hughes  (20)  give  complete  discussions  of  this 
effect.     It  has  not  been  considered  in  this  work. 

Obliquity  of  Foil  with  Respect  to  the  Source 

Even  though  the  direction  of  most  of  the  neutrons  through  the  foil 
does  not  affect  the  activation  directly,    it  does  affect  the  correction 
factors  that  depend  on  foil  thickness  since  the  effective  foil  thickness 
is  a  function  of  obliquity,    (20). 

All  experimental  work  in  this  paper  was  performed  with  the  normal 
to  the  foil  surface  area  facing  the  source.     If  need  arises  to  irradiate 
the  foil  in  a  different  manner,   it  might  be  of  interest  to  experimentally 
determine  the  magnitude  of  this  effect. 

No  effort  to  further  systematize  these  foil  correction  factors  has 
been  made  in  this  work.     It  might  be  desirable  to  have  a  more  com- 
plete study  of  these  factors.     A  determination  of  the  absolute,   unper- 
turbed flux  from  absolute  foil  counting  would  require  application  of 
many  of  the  factors  introduced. 
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The  Kansas  State  University  Graphite  Pile  was  calibrated  for 
absolute  thermal  neutron  flux  with  a  one  curie  Pu-Be  neutron  source, 
thus  standardizing  indium  foils,     a  BF3   probe  and  a  scintillation  probe. 
The  sensitivities  of  these  detectors  to  thermal  and  indium  resonance 
neutrons  are  given. 

Based  on  the  neutron  emission  rate  of  the  source,    absolute  values 
of  the  measured  slowing  down  density  at  indium  resonance  energy 
were  obtained  from  the  volume  integral  of  the  slowing  down  density 
over  a  pile  of  assumed  infinite  height.     The  absolute  thermal  neutron 
flux  was  then  obtained  from  the  thermal  diffusion  equation  using  as  a 
source  condition  the  absolute  slowing  down  density  at  indium  resonance 
energy  extended  to  the  absolute  slowing  down  density  at  thermal  energy 
by  use  of  Fermi  slowing  down  theory. 

The  standardized  detectors  were  used  to  test  experimentally  the 
validity  of  the  diffusion  equation  in  a  finite  pile  geometry,    using  as  a 
source  condition  to  the  thermal  diffusion  equation:     (a)  the  solution 
to  the  Fermi  age  equation  evaluated  at  thermal  energy,     (b)  a  thermal 
point  neutron  source,    and  (c)  the  analytical  expression  for  the  slowing 
down  density  at  thermal  energy  as  found  from  the  standard  pile 
measurements. 

The  techniques  used  in  accurate  foil  activation  measurements  are 
discussed.     Inconsistent  values  of  various  diffusion  parameters  for 
graphite  are  pointed  out.     IBM-650  programs  are  given  to  fit  empiri- 
cally the  indium  resonance  foil  activations,   to  calculate  the  thermal 
flux  from  a  fast  point  source  in  an  infinite  column,   to  calculate  the 
thermal  flux  from  a  fast  point  source  in  a  finite  medium,    and  to  cal- 
culate the  thermal  flux  from  a  thermal  point  source  in  a  finite  medium. 


